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Kurzfassung 


Die ganzheitliche virtuelle Auslegung von Faserverbundbauteilen aus Sheet 
Molding Compund (SMC) ermöglicht es durch die Berücksichtigung von Ferti- 
gungseffekten in der Struktursimulation, bedarfsgerechtere und somit leich- 
tere und günstigere Bauteile zu fertigen. Hierfür muss jedoch die Prozesssi- 
mulation in der Lage sein, die SMC-spezifischen Blockströmung auf Grund der 
niedrig-viskosen Randschicht richtig zu beschreiben. Nur so kann das richtige 
Füllverhalten und somit auch die richtige Faserorientierungsverteilung vor- 
hergesagt werden. Dies wird gerade im Kontext des zunehmenden Einsatzes 
von semi-strukturellen SMC-Formulierungen besonders wichtig. Darüber hin- 
aus gibt es neue Konzepte der Hybridisierung von SMC mit lokalen unidirekti- 
onalen Verstärkungen, welche ebenfalls in der Prozesssimulation berücksich- 
tigt werden müssen, um die Positionierung der Verstärkungen abzusichern. 


Aus diesem Grund wurde in dieser Arbeit ein neuer dreidimensionaler Ansatz 
für die Prozesssimulation von SMC entwickelt, der sowohl die durch Dehnvis- 
kosität dominierte Kernschicht, als auch die niedrig-viskose Randschicht be- 
rücksichtigt. Durch die Verwendung des gekoppelten Euler-Lagrange-Ansat- 
zes, kann auch die Interaktion zwischen SMC und lokalen Verstärkungen 
abgebildet und vorhergesagt werden. Um die Informationen aus der Prozess- 
simulation in die Struktursimulation zu übertragen, wurde eine CAE-Kette 
weiterentwickelt, die das Übertragen der Information vom Prozess- zum 
Strukturnetz ermöglicht. Hierbei wird sowohl eine Cluster-Bildung, als auch 
eine Homogenisierung auf Basis der Faserorientierungsverteilung durch- 
geführt. 


Zur experimentellen Analyse des Fließverhaltens, sowie zur Bereitstellung der 
benötigten validen Materialparameter für die Prozesssimulation, wurde ein 
neues rheologisches Fließpresswerkzeug entwickelt. Mit diesem Werkzeug 
wurden erstmal das kompressible Verhalten von semi-strukturellen SMC For- 


Kurzfassung 


mulierungen nachgewiesen und über einen phänomenologischen Ansatz be- 
schrieben. Auf Basis dieser Kompressibilitätsformulierung konnte ein neues 
rheologisches Model entwickelt, das eine genauere Charakterisierung ermög- 
licht. Durch die Charakterisierung von fünf unterschiedlichen SMC-Material- 
formulierungen, konnten einige Gesetzmäßigkeiten der Materialparameter 
abgeleitet werden. 


Abstract 


A holistic virtual design of Sheet Molding Compound (SMC) components takes 
manufacturing effects in the structure simulation into account. This enables 
lighter and cheaper components that meet the needs ofthe customer due to 
a reduced margin of safety, compared to today’s designs. For that purpose, 
however, the process simulation must be able to describe the SMC-specific 
plug flow correctly, which occurs due to the low-viscosity boundary layer. 
Otherwise, the filling behavior and thus the fiber orientation distribution can- 
not be predicted correctly. Especially in the context of the increasing use of 
semi-structural SMC formulations, this is particularly important. In addition, 
there are new concepts for the hybridization of SMC with local unidirectional 
reinforcements. These reinforcements must also be considered in the process 
simulation to ensure the positioning of the reinforcements in the final part. 


For this reason, a new three-dimensional approach for the process simulation 
of SMC is developed in this thesis. This approach takes into account both, the 
core layer that is dominated by the extensional viscosity and the thin lubrica- 
tion layer, which is dominated by a high shear rate. By using the Coupled- 
Eulerian-Lagrangian (CEL) approach, the interaction between SMC and local 
reinforcements can be modeled and thus, predicted. In order to transfer the 
information from the process simulation to the structure simulation, a CAE 
chain is further developed. This CAE chain enables the transfer of the infor- 
mation from the process to the structure mesh. As part of this CAE-chain, a 
clustering and a homogenization are carried out based on the fiber orienta- 
tion distribution. 


To analyze the flow behavior experimentally and to provide the required valid 
material parameters for process simulation, a new rheological tool is devel- 
oped. With this rheological tool, the compressible behavior of semi-structural 
SMC formulations is proven for the first time and described by using a phe- 
nomenological approach. Based on this compressibility formulation, a new 
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rheological model is developed which allows a more precise characterization. 
By characterizing five different SMC material formulations, several pattern of 
material parameters are derived. 
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1. Preamble 


1.1. Motivation 


In the automotive and in the commercial vehicle industry, long fiber rein- 
forced composites are challenging in terms of two main tasks: weight and cost 
reduction of complex shaped structures. In this context Sheet Molding Com- 
pound (SMC) has an outstanding price to stiffness ratio and a good stiffness 
to weight ratio [Sch07], [Nic11], [AVK13] for medium to high volume produc- 
tions [Äke00], [Sch07]. Due to its possible Class-A surface and its high freedom 
of design, SMC is widely used for exterior parts. By combining SMC with local 
unidirectional reinforcements, the resulting hybrid structure can carry high 
loads in the known load direction, while it can integrate functions and in- 
crease the bending stiffness by adding ribs. 


To design such a complex hybrid structure, a combination of different simu- 
lation steps is necessary. In such a CAE-Chain, the design, the process and the 
structural simulation are combined to consider previous results in further 
simulation steps [Kar15], [Hoh16b]. Since many structural results are strongly 
dependent on the fiber orientation distribution, the process simulation with 
its fiber orientation prediction is an important step in such a CAE-Chain. 


Comparing the results from commercial process simulation tools with exper- 
imental results of compression molding processes, differences in accuracy 
can be noticed, which depend on the polymer type. For thermoplastic mate- 
rials, material models for the injection molding process have been develop- 
ment over a long time. In addition, compression molding simulation of rein- 
forced thermoplastics also shows good agreement with experiments and the 
flow fronts can be predicted even for complex rib structures as shown by 
Hohberg and Baumgartner [Hoh16a]. As the thermoset market is much 
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smaller than thermoplastic injection molding, the prediction of the flow in- 
duced pressure shows significant differences with experimental data due to 
the lack of validated material models [Hoh15]. Therefore, material models, 
which can describe the characteristic SMC plug-flow, are crucial. This plug- 
flow occurs due to the resin-rich lubrication layer with a low viscosity at the 
mold surface. Due to this lubrication layer, no velocity gradient and therefore, 
no stress gradient occurs within the plug flow region. Such plug-flow models 
were developed for 2.5D approaches by Osswald [Oss87], but not for the 
nowadays regular 3-D simulation models. Therefore, no rib structures or T- 
shaped parts can be simulated satisfactorily. Furthermore, the consideration 
of local reinforcements in the molding process needs to be implemented in 
such a new approach, as there is a strong interaction in in-plane as well as in 
the out-of-plane directions. Since these reinforcements can move within the 
tool, the basic assumption for a 2-D or 2.5D approach is not valid any more, 
as velocity gradients through the thickness can occur. 


Material flow properties are important for the correct prediction of material 
flow and the fiber orientation. Since the fiber length of 1 inch (25 mm) is in 
the same range as the typical sample size of shear- or rotational rheometers, 
new experimental rheological measurement approaches are necessary for 
SMC. To consider the material specific flow behavior during the rheological 
characterization, a novel in-line rheometer was presented by Dumont et al. 
[DumO3]. Due to the setup, two separate molding trials had to be performed, 
which caused inaccuracies due to material variations and slightly changed 
boundary conditions like thermal conditions and closing speed and material 
placement. Furthermore, such a setup cannot deal with the compressibility 
of the SMC during the molding trials. To avoid these inaccuracies, an en- 
hanced rheological measurement tool is needed. 


1.2. Objectives 


1.2. Objectives 


The target of this thesis is to investigate and to model the compression mold- 
ing process of Sheet Molding Compound (SMC) with local reinforcements. To 
reach this target, an improvement of the compression molding process and 
especially the rheological material behavior is necessary. Since new SMC ma- 
terial formulations were developed during the last years [Hen17], different 
formulations need to be investigated. Therefore, a new method to character- 
ize the rheological properties has to be developed. This method needs to be 
close to the industrial process to ensure the validity for the further use of the 
rheological parameter. Therefore, industrial standards like the process set- 
tings need to be considered. 


Based on the experimental investigation and characterization, a new SMC 
specific simulation approach will be developed. This approach has to consider 
the two dominating effects of the resin-rich lubrication layer with its high 
shear gradient and the core region with its plug-flow behavior. Since SMC is a 
non-Newtonian fluid, modeling of the shear and extensional viscosity are two 
of the most important steps towards better process simulation. 


To show the capability of this new process simulation approach, a validation 
with a reference structure is performed. Local unidirectional reinforcements 
are also considered in a first approximate model. The results are compared 
with the actual produced part. Furthermore, the fiber orientation results are 
needed in the structural simulation. Therefore, a CAE-chain is developed, 
which can process the result information. This simulation processing output 
contains the export, the mapping and the homogenization, so that structural 
simulation can be performed. 


Briefly summarized, the following sub-objectives can be defined: 


e Experimental investigation of the rheological behavior of different SMC 
formulations 
e Development of a method to measure the rheological properties of SMC 
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e Development, implementation and application of advanced compression 
molding simulation models for SMC with simple local reinforcements 

e Determination of material parameters of different materials for the im- 
plemented models 

e Validation of the implemented models on a reference structure 

e Development of a long fiber CAE-chain to transfer result information 
from the process to the structural simulation 


1.3. Structure of the thesis 


This thesis is divided into four chapters, starting with the presentation of the 
state of research in the area of process characterization and simulation of 
SMC. Therefore, the SMC production and molding process is described in Sec- 
tion 2.1, followed by the modelling approaches in Section 1.2. In Section 2.3, 
the different types of rheological characterization methods for the SMC are 
presented. For the holistic simulation description, fiber orientation develop- 
ment during the molding process and the corresponding closure approxima- 
tions are introduced in section 2.4. In the last section of the first chapter, the 
Coupled-Eulerian-Lagrangian approach is introduced, which enables a fluid- 
structure-interaction in a closed simulation. This is essential for the SMC-tool 
interaction as well as for the SMC-reinforcement interaction. 


The second chapter focuses on the development of a new rheological tool, of 
a new rheological model and on the characterization and evaluation of differ- 
ent SMC material formulations. In Section 2.1, the new rheological tool and 
the experimental setup are presented. In section 3.2, five different SMC ma- 
terials are presented and their experimental characterization tests are per- 
formed. The chosen materials range from the typical Class-A SMC, through 
state of research low-density formulations to three novel semi-structural SMC 
formulations. Based on the experiments, the material compressibility is 
proven for the semi-structural SMCs in section 3.3 and considered for the first 
time in a compressible rheological model in section 3.4. By using this new 
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compressible model, the rheological properties are evaluated in section 3.5. 
By evaluating these properties, physical laws and corresponding modelling 
recommendations could be deduced, which is given in section 3.6. Due to the 
different behavior of the carbon fiber semi-structural formulation, an alter- 
native flow kinetic is presented for such a material in section 3.7. 


In the third chapter, the process simulation of hybrid SMC is presented. As 
basis for development, the state of art software is tested in section 4.1. Inthe 
section 3.2, the new SMC specific 3D simulation approach is developed and 
implemented in the Coupled-Eulerian-Lagrangian method. This is done by di- 
viding the flow into the extensional viscosity dominated core region and the 
shear rate dominated lubrication layer at the tool surface. In section 4.3, this 
model is validated by means of characterization experiments and by a hybrid 
demonstrator, consisting of SMC and local unidirectional reinforcements. 


In the context of the holistic virtual part design, the fourth chapter addresses 
the CAE chain. The export of the process simulation results (e.g. fiber orien- 
tation distribution) is presented in section 5.1. In section 5.2, the mapping and 
in section 5.3, the linear elastic homogenization within the CAE chain is de- 
scribed. At the end of this chapter, an example is given in Section 5.4, showing 
the effects of the consideration of process results in structural simulation. 


2. State of research 


At the beginning of this chapter, the production of the semi-finished Sheet 
Molding Compound (SMC) and the compression molding process is intro- 
duced. This introduction is followed by the state of research for the process 
simulation (Section 2.2) and the rheological characterization methods (Sec- 
tion 2.3) of SMC. For the simulation methods, the state of the art fiber orien- 
tation models (Section 2.4) and approaches for the simulation of material 
flow by using Finite-Element-Methods for large deformations, respectively 
material flow (Section 2.5), are described. 


2.1. Production and molding of SMC 


Sheet Molding Compound (SMC) denotes both, a compression molding pro- 
cess of a stack of planar semi-finished discontinuous fibers and the thermo- 
set-based material itself. The main application for SMC is the exterior of au- 
tomobile and commercial vehicles. Therefore, the semi-finished material is 
mainly based on an unsaturated polyester (UP) resin, which is mixed with 
mineral fillers and other process simplifying and material improving ingredi- 
ents (see Figure 2-1) depending on the designated application [Che16]. Be- 
sides the UP resin, alternative resins like Vinylester resin with higher mechan- 
ical properties, phenol formaldehyde resins and epoxy resins are used 
[Org12], [Hen17]. For a Class-A SMC, the fillers have the highest content in 
the paste. The main task of the fillers is to increase the viscosity and therefore 
to avoid fiber-matrix separation and to reduce the price of the material. De- 
pending on the filler type, it can also functionalize the SMC. For example, Alu- 
minum hydroxide (Al (OH)3) has flame retardant properties, which can lead to 
aviation certified SMCs [AVK13]. To achieve the high Class-A surface quality, 
dispersed thermoplastic additives are added. Due to these additives, the mac- 
roscopic process-induced shrinkage is distributed to many homogeneously 
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distributed micro-areas within the material. Within these micro-areas, the 
thermoplastic additives volume is increasing due to the higher temperature, 
while the curing resin is shrinking. This leads to components which are almost 
completely free of distortion or residual stresses [Kec16]. An important ingre- 
dient is also the thickening agent. It is used to increase the viscosity of the 
paste to handle the semi-finished products and to prevent fiber matrix sepa- 
ration during the molding process [Kec16]. 


Process 
additives 
1,5% 


Figure 2-1: Composition of a typical Class-A SMC in wt.-% [Hoh15] 


All liquid ingredients together with the fillers are mixed and filled in the so- 
called doctor box or knife-coating box (see Figure 2-2) of an SMC conveyor 
belt system. Two boxes are present on an SMC production line to cover and 
coat the upper and lower carrier foils, which are continually moved, with a 
thin resin film. In the next step the continuous fibers, which are fed into the 
chopping units, are cut to the designated fiber length of commonly 1 inch 
(approx. 25 mm). These chopped fibers fall onto the lower foil, covered with 
the resin filler paste, with a random in-plane orientation. By changing the 
speed ofthe foils, the fiber content can be adjusted. After covering the fibers 
with the upper impregnated foil, the whole material is impregnated in a cal- 
endaring zone. In this area, the air is moved out to prevent air entrapments 
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in the final part. In the end of the manufacturing process, the already thick- 
ened material is coiled up or cut into flat sheets [Org12], [AVK13], [Che16]. 


For further detailed information regarding the SMC ingredients and a detailed 
description of the manufacturing process of the semi-finished SMC, please 
see the chapter “Sheet Molding Compound” from Nicolais and Barzachiello 
[Nic12]. 


Top release film Knife-coating box with resin/ Glass rovings 


curing agent/filler mixture 


B 


Roving chopper 


Glass fibers 


Fulling/impregnation belt 
a Bottom release film 


Knife-coating box with resin/ 
SMC take-up curing agent/filler mixture 


Figure 2-2: Principle of SMC production [Che16] 


After a maturing time between two and four weeks, in which the viscosity of 
the semi-finished product is increasing, the material is ready to be processed 
inthe compression molding process. Due to the ongoing increase of viscosity, 
the material is usable between a few weeks up to several month depending 
on the storage conditions, the resin type and the amount of thickening agent 
[Che16], [Org12]. 


For the molding process, the coiled SMC semi-finished material is cut and the 
single sheets are stacked into an initial charge (IC). The number of sheet layers 
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as well as the shape and size of the initial charge depend on the mold size and 
shape (see Figure 2-3). Typically, the initial charge covers between 40 and 70% 
of the projected mold surface, which is also the ratio of the initial filled to the 
total surface are. This reduces the press pressure on the one side, but on the 
other side, also to generate some material flow to release entrapped air and 
process induced styrene bubbles. The molding tool is tempered between 140 
and 160°C depending on the resin system, while the initial charge is placed at 
ambient temperature conditions into the mold. After placing the initial charge 
into the mold, the mold is closing with a closing speed between 1 and 10 
mm st. During the compression molding process, the pressure is increasing 
up to 200 bar for glass fiber SMC and up to 300 bar for carbon fiber SMC. 
These limits are ensured by limiting the maximum press pressure. If this limit 
is reached, the closing speed is reduced and therefore the pressure within the 
mold. After the mold is filled, the material stays in the mold under a high pres- 
sure until most of the curing process is finished. The curing time within the 
mold is between 30 seconds and 3 minutes depending on the resin system, 
the mold temperature and the part thickness. After the curing phase, the final 
part can be ejected [AVK13], [Che16], [Org12]. 
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Parts placed in the mold 


Bottom mold 


SMC components 
Stacked SMC cuttings 


SMC prepreg 


Figure 2-3: Process steps for the manufacture of SMC components with cutting, stacking, 
placing of the initial charge in the mold and the final part [Che16] 


For further detailed information please see the chapter “Sheet Molding Com- 
pound” from Nicolais and Barzachiello [Nic12]. 


2.2. Modelling of the SMC flow process 


For all modelling approaches of the SMC flow process, it can be assumed that 
no curing reaction occurs during the compression molding as long as the SMC 
material is flowing. This assumption is supported by many experimental re- 
sults from Marker and Ford [Mar77], as well as Barone and Caulk [Bar79]. Also 
newer fast curing SMC formulations fulfill this assumptions like Bücheler 
[Büc17] showed in his PhD thesis. Therefore, it is sufficient to describe the 
flow phase by using (thermo-) mechanical models. These models can gener- 
ally be divided into two types with regard to their level of detail. On the one 
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side, there are specialized models, which simulate the fibers and the paste 
separately, so that fiber matrix separation effects (changing of local fiber vol- 
ume content) can be described. These fiber matrix separation effects were 
observed in T-shaped molds [Kim92], [Kim97], [Kuh18] and especially at ribs 
[Xu93], [Chr97]. To model these effects, the fibers are resolved individually to 
simulate the fiber-fiber and the fiber-fluid interactions [Pér15], [Kuh18]. Due 
to the complexity of these models, only small spatial problems can be simu- 
lated with such an approach. Another, less computationally intensive, ap- 
proach is a two-phase model based on the mixture theory on a macroscopic 
scale. This approach describes the two phases, fiber and paste (resin and 
filler), by a volume fraction partitioning. To model the flow of both phases, 
additional balance and constitutive equations are required. Therefore, a mod- 
ified Darcy’s law was used by Orgéas et al. [Org06], [Org07] and Dumont et al. 
[Dum12]. Up to now, these models are just empirical modifications of Darcy’s 
Law, where numerous non-physical model parameters need to be fitted to 
experimental results with high computational effort. 


On the other side, macroscopic one-phase models are used more widely. The 
SMC material is described as a continuum fluid so that the matrix and the 
fibers have the same velocity and therefore no fiber-matrix separation will 
occur. The fiber orientation is then computed from fiber orientation models 
based on physical and empirical models and the flow field (cf. section 2.4). 
The first one-phase models were taking advantage of the shell-like structure 
of the most SMC parts. Hence, the thickness is relatively small compared to 
the other in-plane dimensions. Therefore, generalized Hele-Shaw shell mo- 
dels were developed [Sil80], [Hie80], [Lee84], [Jac86], [Oss90], [Adv90b], 
[Twu93]. The generalized Hele-Shaw shell model describes a 2D flow in a do- 
main with a small thickness compared with the in-plane dimensions. These 
models employ a Newtonian description of the SMC material and the no-slip 
condition between the SMC and the mold surfaces. This leads to a high shear- 
ing of the fluid across the thickness of the mold cavity. 
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Barone and Caulk [Bar85] showed with their experiments the need for better 
and suitable material models. By using different colored SMC sheets and per- 
forming experiments over a wide range of closing speeds, they showed a plug- 
flow behavior of the SMC core region and the presence of a resin rich lubrica- 
tion layer (see Figure 2-4). Furthermore Lee and Tucker [Lee87] showed the 
dependency of the lubrication layer thickness on the compression molding 
process settings. 


Figure 2-4: Experimental investigation of the flow behavior of SMC by using differently col- 
ored SMC sheets [Bar85] 


Based on these experiments considering the lubrication layer, analogue ma- 
terial models were developed. In a first step, different friction formulations 
were investigated to consider the lubrication layer in a generalized Hele-Shaw 
shell model by Barone and Caulk [Bar86]. By comparing the friction mecha- 
nisms of constant friction, hydrodynamic friction and Coulomb friction in the 
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simulations with the experiments, the hydrodynamic friction approach based 
on the relative velocity at the mold surface gave the best consistency. Based 
on this result and on the consideration of uniform in-plane shear and exten- 
sional deformation [Kot98], [DumO3] as well as only out-of-plane shear in the 
lubrication layer [Ode04], various new material flow models were developed 
and improved [Bar87], [Bar88], [Cas89], [Oss94], [Abr03], [Dum03], [Dum07], 
[Gui10], [Gui12]. These models use the superposition of the fluid pressure p, 
the viscous stress of the core region Qype. and the friction stress from the 
lubrication layer Of,i¢, by assuming an incompressibility of the material: 


= —pl + Orheo + Ofrics (1) 


were the incompressibility for a control volume is defined as: 


Vy Ivy Ovz _ 

Oe Dr + DE = 0 and (2) 
d 

= p(t x(t) = 0, H 


were 1; is the velocity in x, y or z, and p the density of the material. 


As a result of these one-phase models, which take into account the core layer 
and the lubrication layer behavior, it can be stated that certain modelling ap- 
proaches have proven their effective value. As the lubrication layer is rela- 
tively thin compared to the part thickness, the hydrodynamic friction ap- 
proach by using a power-law model becomes the state of the art in research 
studies [Abr03], [Dum07], [Org12], [Gui12]: 


Ofric = —A Op = Qn. (4) 


vo vll 


In this equation, represents the velocity in the surface plane, A is the hydro- 
dynamic friction coefficient, m is the power-law coefficient and N the normal 
surface vector. Furthermore, ® represents the dyadic product of both vec- 
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tors. In this equation, the velocity ¥ and the surface normal vector N are per- 
pendicular to each other. Therefore, the symmetry of the friction stress ten- 
SOr Orr. is scarified. 


For the viscous stress of the core region, a deformation rate dependent 
power-law approach becomes the state of art in research [Le 02], [Dum03], 
[Dum07], [Org12], [Gui12]: 


IDIN” 

Cno = 4 (5) D- 5) 
0 

In this equation, 7 represents the extensional viscosity for a pure strain load 

case, D the deformation rate tensor, which is deviatoric, and n the viscous 

stress power-law coefficient. 


Motivated and based on the experiments by Barone and Caulk [Bar85] a mul- 
tilayer plug flow model, considering the single SMC sheets, was developed by 
Leterrier and G’Sell [Let96]. This model is also based on a generalized Hele- 
Shaw approach, but calculates the deformation of every single sheet and their 
interactions. As the result is similar to the previous Hele-Shaw shell models 
but more comprehensive, this approach was not continued. 


None of the presented one-phase approaches considers the anisotropic flow 
properties of the SMC material due to its fiber orientation. Depending on the 
process condition (e.g. initial charge, closing speed), anisotropic material 
properties can be observed in the final part structure and also a change in the 
flow was assumed by Kim et al. [Kim92]. Therefore, a Hill’s anisotropic yield 
criterion in conjunction with a rigid visco-plastic formulation was chosen and 
a resulting finite element formulation was developed [Kim92], [Lin97], 
[Lin99]. These models were able to describe the slight changes due to the an- 
isotropy, at the expense of six additional parameters, which must be deter- 
mined as a function of the strain rate. Due to the complex modelling and char- 
acterization, this approach was not continued. 
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Since all these approaches are based on a 2.5D model, a limitation in the pre- 
diction of complex shaped parts can be identified. For example in ribbed 
parts, the Hele-Shaw boundary condition of a relative low height compared 
with the in-plane 2D dimension is not satisfied anymore. Therefore, a new 3D 
approach is developed in this thesis to overcome the Hele-Shaw limitations. 
By developing such a 3D approach, new functionalities like the interaction 
with local reinforcements or inserts are possible. Therefore, the state of the 
art research approaches from the 2.5D simulations are used and extended to 
the 3D simulation. 


2.3. Rheological characterization of SMC 


The aim of the rheological characterization of SMC is to determine the neces- 
sary material parameters for the flow models during the compression mold- 
ing process. This can be done by two different methods. On the one side, 
shear or rotational rheometers can be used. This is an established method to 
characterize the rheological properties of polymers, but has the disadvantage 
that the induced flow in a rheometer is different to the flow in a mold. Ina 
rheometer, high shear rates are applied, while SMC shows a bulk flow behav- 
ior due to the lubrication layer in the compression molding process (cf. Sec- 
tion 2.2). Furthermore, typical process parameters such as the applied pres- 
sure cannot be realized in a reasonable way. Besides this, the typical sample 
diameter size for a shear or rotational rheometer is as long as the fiber lengths 
of the SMC (1 inch, 25 mm) [Lee81]. 


In contrast to traditional laboratory scale rheology, in-line rheological com- 
pression molding tools and corresponding characterization models were de- 
veloped. These in-line rheological tools use different flow kinetics and there- 
fore, different input types. On the one side, in-situ methods are used which 
utilize in-mold pressure sensors, press force or transparent mold sides to fol- 
low and record the flow fronts. These in-situ flow kinetics include squeeze 


16 


2.3. Rheological characterization of SMC 


flow, simple compression, channel flow and simple shear. In addition, a char- 
acterization based on the flow length depending on the process parameters 
was performed using spiral flow. 


The squeeze flow method [Sil81], [Lee81], [Kot98] is very simple (see Figure 
2-5) and one can observe the material over a large range of strain rates by 
varying the closing speed for radial flow of the sample. By recording the press 
force and the sample thickness by the distance between the mold sides, a 
generalized Hele-Shaw model (cf. section 2.2) can be used to determine the 
material parameters. For these tests, it is crucial that no slip between the SMC 
material and the mold surface occurs, as this would lead to an extra shear 
stress, which would affect the measurements. Due to the setup, the flow kin- 
ematic is neither homogeneous over the sample thickness nor along the flow 
radius. Therefore, the inversion of the flow models from section 2.2 needs to 
be done iteratively to consider the correct flow kinematics while fitting the 
material parameters. This can easily lead to errors in the material parameters. 
Therefore, these types of rheological testing is mainly done for simple rheo- 
logical models which assume incompressibility of the material [Org12]. 


FIXED PLATE 


Figure 2-5: Schematic illustration of the squeeze flow rheometer [Sil81] 
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The simple compression in-line rheometer [Mic90], [Le 02], [Dum03], 
[Dum07], [Mar09], [Gui12] is based on the same idea as the squeeze flow rhe- 
ometer, but overcomes its disadvantage by lubricating the mold surfaces and 
using a higher initial radius to height ratio (see Figure 2-6). Due to the slip 
condition, as a result of the lubricated mold surface, the material flow is more 
homogeneous. As the experiments were performed with constant strain 
rates, the extensional viscosity can be calculated directly by using eq. (5). As 
the lubrication layer is suppressed in these experiments, the friction parame- 
ter cannot be determined. Furthermore, the material flow in the rheometer 
is different to the flow in the real tool, as the friction stress is missing. 


1- Fixed plate 
2- Moving plate 
3 - SMC disk 

4 - Lubricant 


5 - Heating cartridge 


Figure 2-6: Schematic illustration of the simple compression rheometer. To gain simple 
compression, the mold surfaces are lubricated [Dum03] 


The plane strain/channel flow rheometer [Cas90], [Lin97], [Lin99], [Dum03] 
simplifies the flow compared to the simple compression rheometer. Instead 
of molding the sample in a 2D flow, the material flow is reduced to a quasi 1D 
(see Figure 2-7). These tests were performed with and without lubricating the 
mold surfaces. Due to the simplification of the flow and the high ratio of the 
initial charge length to the initial charge height, no flow kinetic gradient oc- 
curs within the material. Therefore, a generalized Hele-Shaw shell model with 
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a deformation rate based core region and the lubrication layer can be devel- 
oped by using the eq. (1) - (5). By inverting this model, the material parame- 
ters for both regions can be determined. 


1- Fixed plate 
2 - Heating cartridge 


3- Moving plate 
4 - SMC sample 


5 - Lubricant 


6 - Transverse forcesensor 


7 -Slide 


he 


Figure 2-7: Schematic illustration of the plane strain/channel flow rheometer. To gain 
plane strain, a 1D flow is enforced [Dum03] 
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The in-line rheometer based on the simple shear kinetics [Le 02] is an up- 
scaled rotational rheometer. Due to the larger scale, a sufficient ratio be- 
tween fiber and sample diameter as well as sample diameter to height is en- 
sured. To get closer to the industrial molding conditions, normal forces are 
applied. This also helps to prevent the slip between the material and the mold 
surface. In some cases, this is not sufficient to prevent slip and an additional 
adhesive film can be introduced between the mold surface and the sample. 
Due to the no-slip condition of the material at the mold surface, these exper- 
iments could only be performed at room temperature to prevent the devel- 
opment of the lubrication layer, as higher temperatures would lead to a slip. 
Due to the plug flow behavior of the core region of the SMC, the shear behav- 
ior is mainly important for the hot resin-rich lubrication layer. Due to these 
problems, simple shear rheometers have not been further investigated. 


The spiral flow [Rab08] is mainly used to compare materials under the same 
process conditions. Therefore, the materials are placed in the center of the 
tools (see Figure 2-8) and the flow length along the spiral is compared. By 
changing the process conditions, the best material for specific molding pro- 
cesses can be identified quantitatively. Besides this method, the flow along 
the spiral can be approximated as a 1D flow. Therefore, the tool can be used 
as a plane strain rheometer. By using the pressure sensors and the press pres- 
sure, the hydrodynamic friction coefficient can be calculated. By comparing 
the results between this spiral flow rheometer and the plane strain rheome- 
ter, a good match could only be achieved for low fiber volume contents and 
medium filler contents. These differences might come from the approxima- 
tion of the plane strain, which is not completely satisfied in this case. There- 
fore, this type of rheometer has not been investigated for the characteriza- 
tion of SMC materials regarding the process simulation. 
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Figure 2-8: Schematic illustration of a spiral flow tool [Rab08] 


The plane strain rheometer becomes state of the art in research for the rhe- 
ological characterization of SMC. They are close to the industrial process and 
therefore, the characterization results can be used for industrial applications. 
The available plane strain rheometers have the huge disadvantage that at 
least two measurements with different initial charge positions are needed. By 
comparing the different pressure distributions of these different initial charge 
positions, the characterization is performed. Due to the reproducibility of 
flows and the local variation of the SMC material properties, this leads easily 
to small inaccuracies due to experimental errors. To overcome this disad- 
vantage, a new plane strain rheometer is developed in this thesis, with which 
a characterization within a molding trial could be performed. Due to the new 
design, the flow front can be traced during the molding process at several 
points. This enables additional characterizations of the material flow and 
therefore, a more precise understanding of the SMC compression molding 
process. 
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2.4. Modelling of the fiber orientation distribution 


As explained in section 2.2, the main macroscopic models describe the SMC 
as a one-phase fluid. Therefore, the fibers and their orientation are modeled 
as a function of the flow. Such a homogenization of particles were firstly in- 
troduced by Einstein [Ein1906]. Based on this idea, Jeffery [Jef22] developed 
the first fiber orientation model for the rotation of a single ellipsoid fiber. By 
assuming a Newtonian fluid, a linear velocity field far away from the fiber and 
no buoyancy or inertia effects, an exact and closed solution has been devel- 
oped. The model describes the two main fiber orientation phenomena: the 
fiber alignment in the shear direction and the alignment in the stretching di- 
rection. The initial orientation of a single fiber at the beginning of a calculation 
step is described by the unity vector p and based on the following expression 
of reorientation: 


2=W-p+&(D-p- [p-D -plp), (6) 


where W represents the vorticity tensor (see eq. 50), D the rate of defor- 
mation tensor and £ the particle shape defined as: 


ar?-1 
ar2+1’ 


g= 


(7) 


where the aspect ratio of the fiber is ar = l/d (& = -1: disc, = 0: sphere, 
€ = 1: ellipsoid with l >> d). 


In the next step, a model to predict the fiber reorientation in concentrated 
suspensions was developed by Folgar and Tucker [Fol84]. Based on their ex- 
periments, they added an empirical rotary diffusion term to the Jeffery’s 
equation (eq. (6)) to model the fiber-fiber interactions. As not all fibers are 
aligned in the same direction, an efficient way to store this information is nec- 
essary. By introducing the fiber orientation function Y(p), which is fitting the 
angular distribution over the in-plan orientation, a new storage method was 
developed: 
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LYP dp = 1, (8) 


where 1. dp denotes the integral over the area of all possible directions of 
the fiber orientation vectors p. To represent a 3D fiber orientation, the fiber 
orientation tensor A was introduced [Ken84], [Adv87]. The fiber orientation 
tensor is a symmetric second order tensor: 


A= [,p®p\W(p) dp, (9) 
or fourth order tensor: 
A= [,P@p@p@p(p) dp, (10) 


with the trace of one, due to the statistical distribution function. An example 
of the fiber orientation distribution and the corresponding in-plane second 
order fiber orientation tensor is given in Figure 2-9. The model developed by 
Folgar and Tucker gave good agreement with short fiber suspensions, but it 
raised the problem with the need of the fourth order tensor (A) to calculate 


i : : d ; 
the second order fiber orientation tensor change A = A). Therefore, a clos- 


ing approximation between the second and the fourth order tensor is needed 
and discussed later in this Section. 
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Figure 2-9: Fiber orientation distribution and the corresponding reduced 2D second order 


fiber orientation tensor 


More precise measurements, considering the orientation distribution over 
the thickness and over the time, showed a development of different fiber ori- 
entations in the core and the shell layer which could not be predicted cor- 
rectly by the Folgar-Tucker model. Based on the idea of the Strain Reduction 
Factor (SDF) which describes the local reduction of the strain due to resin rich 
“slip layers” with higher strains, the Reduced Strain Closure (RSC) was devel- 
oped by Wang et al. [Wan08]. Therefore, the scalar reduction factor k was 
introduced, which has a value between zero and one and needs to be fitted 
to experiments. 


At the same time, new developments were done for the behavior of long fi- 
bers. In the basic Folgar-Tucker model, the interaction coefficient is a scalar 
value that models isotropic rotary diffusion. This is suitable for short fibers, 
but as long fibers cannot be assumed to be rigid, the interaction depends on 
the current fiber orientation. Therefore, the scalar value is replaced by a ten- 
sor, which can vary between the interaction of aligned and quasi-isotropic fi- 
ber orientations. Furthermore, this Anisotropic Rotary Diffusion (ARD) tensor 
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is flow type dependent (shear and extensional) [Phe09a]. This tensor contains 
five new parameters (bı - bs), which need to be fitted to experiments. 


As the RSC and the ARD extension is dealing with two separate effects of the 
fiber orientation evolution, the combination of both should be used for long 
fiber reinforced polymers such as SMC. The combined ARD-RSC Folgar-Tucker 
model has the form: 


A=(W-A-A:-W)+E(D-A+A-D-2[A+ (1- K)(L-M: A)]: 
D) +y(2[C — (1- K)M:C]—2xtr(C)A—5(C-A+A-C) + 10[A+ 
(1- x)(L—- M: A)l:C), (11) 


with the ARD fiber interaction tensor: 
_ 2 D D? 
C= b,I+ b,A+ b,A +b,+ bs 7z (12) 


and the fourth order tensors L and M, which are analytical functions of the 
eigenvalues and eigenvectors of the second order fiber orientation tensor 
A. They are defined as 


L= 1 Ale; Qei © e; Q e;) (13) 
and 
M= rhe; Qe We; 8 e;). (14) 


an 


Furthermore, the “:” represents the double dot product of the second and 
forth order tensors. 


For the closure from the second fiber orientation tensor A to the fourth order 
tensor A, different approaches have been developed and tested. The linear 
closure approximation [Han62] gives exactly the fourth order tensor for an 
isotropic fiber orientation, but the accuracy is getting worse when it comes to 
aligned fiber orientations. The quadratic closure approximation developed by 
Doi [Doi81] and Lipscomb et al. [Lip88] uses the dyadic product of the second 
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order tensor with itself. This gives a good agreement for aligned fibers, but 
shows some nonphysical oscillation effect depending on the flow regime. 
Therefore, the linear and the quadratic closure approximation were com- 
bined in the hybrid closure approximation by Advani and Tucker [Adv90a]. 
Depending on the fiber orientation, a scalar value is used to interpolate be- 
tween both approaches. 


By showing the existence of an exact closure, if fiber-fiber-interactions are 
neglected [Lip88], the natural closure approximation was developed by Ver- 
leye et al. [Ver94]. This method uses a general expression for the fourth order 
tensor, which can be calculated by an analytical solution. Therefore, the nat- 
ural closure values need to be fitted. This method showed a good agreement, 
but has some problems with singularities. 


The orthotropic fitted closure approximation was developed by Cintra and 
Tucker [Cin95] and uses the assumption that the second and the fourth order 
tensor have the same principal direction and that the principal values of the 
fourth order tensor are a function of the second order tensor’s eigenvalues. 
The parameters for this correlation were also fitted from experimental data 
and obtained from distribution function calculations. Due to nonphysical os- 
cillations an improved version (ORW3) was proposed by Chung and Kwon 
[Chu01]. Due to the high number of parameters and the eigenvalue calcula- 
tions, this method is very computational intensive. 


To find a closure approximation which combines the accuracy of the eigen- 
values methods with the relative low computational time of the natural clo- 
sure approach, the Invariant Based Optimal Fitting closure approximation 
(IBOF) was developed by Chung and Kwon [Chu02]. This method uses a poly- 
nomial function of the invariants of the second order fiber orientation tensor. 
Therefore, the components of the fourth order tensor can be calculated us- 
ing: 


Aij = Bi S(5i; Ox) + Bo S(6i; Au) + Bs S(Aij Arı) + 
Ba S(6i Akm Amt) + Bs S(Aij Akm Amt) T Bo S(Aim Anj Akn An) , (15) 
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where S is a distribution function of the components in the fourth order Ten- 
sor A, B3, B4 and Be are coefficient depending on the second and third invari- 
ant of the second order fiber orientation tensor and ßı, B2 and Bs are coeffi- 
cients depending on the invariants and B3, B4 and Be. For further details on 
these coefficients, it is referred to Chung et al. [Chu02]. As the coefficients B3, 
Ba and Be are polynomial functions, each coefficient has 21 parameters. A to- 
tal number of 63 parameters need to be fitted based on six known exact so- 
lutions of the fourth order tensor. This fitting was done for different thermo- 
plastic short and long fiber materials by Phelps [PheO9b]. The parameter set 
for the long fiber reinforced thermoplastics (LFT) is used for SMC as well, since 
the fiber length and fiber volume fractions are in the same range. 


Summarizing, the ARD-RSC becomes state of the art for long fiber reinforced 
polymers, since it has the best agreement with experiments. Furthermore, it 
has an acceptable parameter identification effort by fitting the five coeffi- 
cients of the ARD fiber interaction tensor C and the scalar reduction factor k. 
For the closure approximation of the fourth order fiber orientation tensor A, 
the IBOF closure is a good compromise between accuracy and computational 
time. The parameter determination for the IBOF closure is quiet complex, but 
due to the available parameter set, the combination of the ARD-RSC model 
with the IBOF closure is used in this work. Since the fiber orientation model- 
ling is not the focus of this thesis, no further developments are done in this 
context. 


2.5. Coupled-Eulerian-Lagrangian approach 


To model a flow of material, a Lagrangian mesh based method can be used 
with some difficulty. As a high deformation is necessary, Lagrangian meshes 
might be distorted during the deformation and the Jacobian determinants 
might become negative. This could be overcome by remeshing the Lagrangian 
mesh, but this approach leads to a higher computational time, can introduce 
errors from the projection and can cause convergence problems. Therefore, 
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mainly Eulerian mesh formulations are used for fluid transport problems. Due 
to the convection of mass through the elements, the constitutive equation 
and its update can become complicated. Furthermore, the consideration of 
moving boundaries and interfaces is difficult. To overcome these drawbacks, 
a hybrid method was developed which combines the best of the Eulerian and 
the Lagrangian methods while minimizing their disadvantages. The Arbitrary 
Lagrangian Eulerian (ALE) method was first introduced by Noh [Noh63] and 
improved by Hirt et al. [Hir74]. Belytschko and Kennedy [Bel78a] showed the 
advantages of these ALE methods for three different use cases in 2D and 3D. 
They also showed the capability of the ALE formulation in terms of the fluid- 
structure-interaction [Bel78b] and Hughes et al. [Hug81]. This ALE method is 
distinguished by the mesh motion, which is separated from the material mo- 
tion. Due to this separation of the material and the mesh motion, a material 
flow through the ALE mesh as well as a mesh deformation can occur. Depend- 
ing on this mesh motion , the ALE approach is similar to the Eulerian method, 
if the mesh motion is set to zero, while it describes a Lagrangian formulation, 
if the mesh motion is equal the material motion [Bel14]. Furthermore, differ- 
ent dynamic, stress dependent mesh motions were investigated by Liu et al. 
[Liu91]. 


In this general ALE method, the material domain, the reference domain 
(mesh) and the spatial domain have to be described separately by defining 
the motion of the materials as: 


x= P(X,t), (16) 


with the material coordinate X. Therefore, the function @(X,t) maps the 
body from the initial configuration Qo to the spatial domain Q. Since both do- 
mains are not directly linked to the mesh deformation, the reference domain 
Q is introduced in the ALE Continuum Mechanics (see Figure 2-10). Within this 
domain, the initial values of the particle motions X are denoted as: 


X= $(X, 0). (17) 
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By introducing this reference domain Q), the motion of the material can be 
described mesh-independently. In most implementations, this domain is used 
to construct the initial mesh. Since the reference domain Î typically remains 
coincident with the mesh, it is also the computational domain. To map the 
material motion from this reference domain to the spatial domain Q, the 
function 


x= dw) (18) 


is used. Therefore, the kinetics for the general ALE formulation and their sim- 
plifications towards the Lagrangian, and the Eulerian formation is given in Fig- 
ure 2-11. For further information regarding the general ALE formulation, see 
the chapter “Arbitration Lagrangian Eulerian Formulations” in [Bel14]. 


Spatial Domain N 
.x 


Reference Domain 0. Material Domain No 
ae -x 


Figure 2-10: Maps between Lagrangian, Eulerian and ALE domain 
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Description General ALE Lagrangian Eulerian 
Motion Material x=0(X, 1) x = 0(X, f) x = O(X, t) 
Mesh x= O(X, 0) x= (X, 1) x= I(x) 
(X=X,6=9) (X =x, =I) 
Displacement Material u=x-X u=x-X u=x-X 
Mesh ü=sx-y ü=sx-X=u ü=x-x=0 
Velocity Material V=U, iy VEU, iy) V=U, m 
Mesh v=ü,, t= ûn ¥=6,,.,=0 
Acceleration Material a=v, in a=v, ayy a=v, iy 
Mesh a= Yny a=V,y=2 à=? 
Figure 2-11: Comparison of the kinetic for an ALE formulation with Lagrangian and Eulerian 


description [Bel14] 


The Coupled-Eulerian-Lagrangian (CEL) Method is a special ALE formulation, 
in which the Lagrangian formulation is combined with an ALE formulation (see 
Figure 2-12). This is done to improve the interaction of the flowable material 
with the Lagrangian solid material. Therefore, the contact and moving bound- 
ary conditions are treated in the Lagrangian step. By coupling an ALE step with 
an Eulerian element formulation to this step, high deformations can be real- 
ized. By using a negative mesh motion in this ALE step, which is reverse to the 
material motion, the mesh stays the same. As a consequence, the final con- 
figuration is comparable to the Eulerian formulation. Different methods of 
remeshing were investigated by Huerta and Liu [Hue88]. For this remapping, 
they use an approximation of the fluid boundary surface in the deformed con- 
figuration. Based on this idea, various further works were focusing on the re- 
duction of the computational effort due to new terms in the formulation for 
special geometries [Ben89], [Ben92]. Besides the numerical remeshing, the 
material point information needs to be transported due to the material con- 
vection. Therefore, Liu et al. [Liu86] developed the convection equation to 
fulfill the mass conservation by using the Lax-Wendorf method [Lax60]. 
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For further information regarding the operation of the Eulerian, the Eulerian 
operator splitting and the Lagrangian step see Benson et al. [Ben04] and Do- 
nea et al. [Don04]. 


ee Step Eulerian Step 
ee 


Figure 2-12: Lagrangian and Eulerian Steps in the CEL approach [Ben04] 


The interaction in the Lagrangian step is done by using a penalty function. 
Therefore, the penetration length g between the element edges or the inter- 
polated surface of the material in the Eulerian elements and the interaction 
surface is calculated (see Figure 2-13). The resulting penalty force depends on 
the exact formulation of the penalty function, but it is defined normal to the 
interaction surface towards the material. For an unambiguous description of 
this penalty force direction, it is recommended to use solid interaction ele- 
ments as they have a distinctive surface direction. Furthermore, at least three 
Eulerian elements behind the Lagrangian contact surface are needed to pre- 
vent a loss of mass over insufficient boundary conditions. 
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Rigid Body (tool) 


Rigid Body (tool) 
4 


Pure Element Void Element 


Figure 2-13: Definition of the penetration for pure and void elements [Ben04] 


The current focus of research regarding the ALE and CEL methods has two 
main targets: the simulation of cracks in a structural simulation and the de- 
velopment of special material models in the fluid simulation. Since the ALE 
and CEL formulation are not limited to fluids, research for high deformation 
structural simulation was also done. Currently, combination with XFEM meth- 
ods leads to new application areas. For example the prediction of crack prop- 
agations due to contact is simulated by Vitali and Benson [Vit08]. On the other 
side, the ALE and CEL formulations are improved for special applications and 
material combinations. Therefore, a multi-physics and multi-material code is 
developed to simulation the interaction of plasma to solid materials by Ko- 
niges et al. [Kon15]. 


In this thesis, the commercial available CEL framework within Abaqus is used. 
Therefore, no further developments or implementations of the CEL formula- 
tion need to be done. Since the CEL implementation within Abaqus offers 
most of the user subroutines, the focus is on the material model development 
and the corresponding CEL specific simulation parameter identification. 
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3. Rheological measurements and 
rheological modelling 


In this chapter the rheological measurements, the rheological modelling and 
the evaluation of different SMC materials is presented. Therefore, the meas- 
urement method is described in Section 3.1. In the next Section, the experi- 
mental investigated are presented and measurement results are shown. In 
Section 3.3 and 3.4, the compressibility and rheological models are devel- 
oped. In Section 0, the correlation between the material parameters for the 
different materials and their composition is elaborated. 


Most of the content of this chapter 2 has been published predominantly in 
the publications [Hoh17a], [Hoh17b], [Hoh17c]. 


3.1. Rheological tool and experimental setup 


As discussed in Section 2.3, typical shear or rotational rheometers cannot be 
used for the rheological characterization of SMC, since the fiber length of SMC 
is in the same range as the typical sample size. Furthermore, SMC has a dif- 
ferent flow behavior compared to other thermoplastic and thermoset com- 
posite materials in the injection molding or compression molding process (cf. 
Section 2.2). Therefore, not all necessary material parameters can be deter- 
mined by using a classical rheometer method. 


The first characterization of SMC by using an SMC tool was performed by Cas- 
tro and Griffith [Cas89] and improved by Dumont et al. [Dum03], [Dum07], 
[Gui12] (cf. Section 2.3). Based on their experiments a new rheological tool 
was developed at the Fraunhofer Institute for Chemical Technology (ICT), 
Pfinztal, Germany. The design of the new rheology tool was driven by the two 
targets of having a long 1D flow of the SMC and the ability to measure the 
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pressure distribution over the flow length. Both targets were deduced from 
the goal to observe the changes in the material properties over the flow 
length and over time. To enable a long flow length, the rheology tool has a 
rectangular cavity of 800 mm x 250 mm (see Figure 3-1) and can produce 
plates between 1 mm and 5 mm thickness. By changing the initial charge, 
mold coverage between 20 % and 80 %, different flow-lengths and their char- 
acteristics can be observed. [Hoh17a] 


For a precise rheological measurement, the knowledge of the time- and loca- 
tion-resolved pressure distributions are important. To gain this pressure in- 
formation, seven pressure sensors are placed in the lower mold part along 
the flow path (e1) (cf. Figure 3-1). They are centered vertically to the flow di- 
rection in the symmetry line of the tool (e2 = 0). Figure 3-1 shows the exact 
position of these sensors, which were mounted flush. These sensors are not 
distributed equally as the tool contains a replaceable inlay in the center area 
and therefore, no sensors could be placed at the transition between the tool 
and the inlay. Furthermore, Sensor 1 and 7 are not equally spaced to the tool 
sides, since the short tool side is replaceable and so the distance needs to be 
increased to protect Sensor 1. 


| ' 800mm ; | 


Figure 3-1: Position of the SMC initial charge and the pressure sensors along the lower 
mold part of the rheological tool [Hoh17a] 
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Besides the pressure information, the press data like the current press work- 
ing force and the position of the parallelization-equalizing cylinders are nec- 
essary for the rheological characterization. Therefore, the internal press sys- 
tem is linked to the same data logging system as the pressure sensors. With 
this experimental setup, it is possible to record all necessary data simultane- 
ously (see Figure 3-2). 


This experimental setup was used at two Dieffenbacher presses at the Fraun- 
hofer ICT. The bigger Dieffenbacher COMPRESS PLUS DCP-G 3600/3200 AS 
has a maximum working pressure of 36 000 kN and a maximum closing speed 
of 80 mm s”. The smaller Dieffenbacher DYL 630/500 has a maximum work- 
ing pressure of 6 300 KN and a maximum closing speed of 18 mm s“. 


Figure 3-2: Experimental setup for the rheological measurements 


3.2. Rheological measurements of different SMC 
materials 


During the experiments, materials with different resin types, fillers and fibers 
were investigated. All these SMCs were compounded at Fraunhofer ICT in 
Pfinztal, Germany, on a modified HM-LB-800 Schmidt & Heinzmann SMC line 
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(cf. Section 2.1). The detailed description of the five investigated SMC formu- 
lations (Class-A SMC, the Low-Density Class-A SMC, Vinylester SMC and the 
UPPH carbon and glass fiber SMC) follows within this Section. For each SMC 
formulation a short description of the material composition is given. This is 
followed by a description of the process parameters and a presentation of a 
reference pressure distribution plot for the same closing speed of 


h= 1mms™t. 


3.2.1. Standard Class-A SMC 


This first SMC, which is investigated, is a typical Class-A SMC for exterior parts 
with a paintable surface. It has a high amount of fillers to reduce shrinkage 
and material costs while increasing the density only slightly. The resin-filler 
mixture (see Table 3-1) has a density of pciass-A, paste = 1.69 g cm, The semi- 
finished product was produced with a fiber content of 20 vol.-% chopped glass 
fibers (27 wt. %) with a fiber length of 1 inch (25 mm). The fibers used for the 
SMC are a 12k MultiStar® 272 from Johns Manville (filament diameter 
d = 13.5 um, 4800 tex). As filler, calcium carbonate (CaCOs, chalk) is used. Due 
to its low price, this is a typical filler for Class-A SMC formulations (see Section 
2.1). The final semi-finished material has a density of pciass-a = 1.85 g cm. The 
exact composition of this Class-A SMC is given in Table 3-1. 
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Table 3-1: | Composition of the paste of the Class-A SMC 


Component Trade name Quantity 
UP resin Palapreg Premium G21 100 parts 
Adherent and flow aids BYK 996 + BYK A560 2.2 parts 
Peroxide Pergan Peroxid BEC 1.5 part 

L&V 50%MgO Luvatol® MK35 2.5 parts 


Filler: calcium carbonate 


én r 
(CaCO3) Omya Millicarb® OG 190 parts 


The experiments with this Class-A SMC were performed using an initial mold 
coverage of 37.5 % (300 mm x 250 mm). The room-tempered initial charge 
was assembled out of four single sheets and had an average thickness of 
ho= 8.4 mm. The mold was tempered to 150 °C for the upper side and to 
160 °C for the lower side. Due to the different thermal expansions at the up- 
per and lower tool sides, this helped to prevent both tool sides to touch each 
other at the flush faces. Furthermore, these settings helped to balance tem- 
perature differences due to convection. [Hoh17b] 


All compression molding experiments for this material were performed with 
a constant closing speed of h = 1 mm st, which was predefined directly in the 
press control. To prevent too high pressure in the cavity and therefore, to 
protect the pressure sensors, the maximum working force was limited to 
2 000 kN. This leads to a reduction of the closing speed at the end of the com- 
pression molding process at a relative gap height of 0.25 mm (see Figure 3-3). 


The pressure distribution and the press working force over the relative gap 
height (h - ht) is plotted in Figure 3-3), where h is the current height and hr is 
the final height of the plate. This type of plotting enables an easier compari- 
son of the slightly different flows and initial charges. [Hoh17b] 
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Figure 3-3: Exemplary pressure distribution of the standard Class-A SMC for a closing speed 


ofh=1mms! 


The pressure distributions of all 10 measurements are following the same 
pressure pattern. The Sensors 1 to 3 were covered by the initial charge. There- 
fore, these sensors show an immediate pressure increase, as soon as the up- 
per mold side is touching the initial charge material. The Sensor 4 right behind 
the initial charge (in flow direction) can be used to monitor the start of the 
material flow. By comparing the gap heights of the first pressure increase of 
this reference sensor and the one of Sensor 1, an offset between these gap 
heights can be observed. This offset distance can be attributed to a short com- 
pression or squish [Old04] phase at the beginning of the mold filling process. 
After this transversal compression (e3), the flow of the SMC starts and the 
pressure at Sensor 1 reaches a threshold at a relative gap height of approx. 
5 mm. During the following stable plug-flow, the pressure at Sensor 1, 2 and 
3 is increasing throughout [Hoh17b]. After the press working force reaches its 
maximum, at a relative gap height of approximately 0.25 mm, the pressure 
values of the sensors are decreasing due to the reduced closing speed. 
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3.2.2. Low-Density (LD) Class-A SMC 


The second SMC formulation is a Low-Density (LD) Class-A SMC. This LD Class- 
A SMC has a lower density compared to the previous standard SMC due to 
micro hollow glass spheres, which have a diameter of 40 um and a density of 
0.38 g cm”. The formulation was developed by the Fraunhofer ICT during the 
Technologie-Cluster Composites “InnoSMC” [Kec16]. Using the compo- 
sition from Table 3-2, the density of the paste was reduced to 
PLD-Class-A, paste = 1.09 g cm”. The resin types for the Class-A and the LD Class-A 
SMC are not the same, as the resin producer was not offering the same resin 
type anymore. By evaluating the successor resin, no differences could be 
found and therefore, the resins can be considered to have similar behavior. 
The semi-finished SMC was produced at the Fraunhofer ICT with a fiber con- 
tent of 20 vol.-% (38 wt.-%) and a density of puo-ciass-a = 1.38 g cm’. As for the 
Class-A SMC, chopped 12 k MultiStar® 272 from Johns Manville with a fiber 
length of 1 inch (25 mm) is used for this SMC. 


This LD Class-A SMC has nearly the same mechanical properties as the stand- 
ard Class-A SMC, while the costs increasing only slightly [Kec16]. Since the 
material still has a Class-A quality at reduced density, it becomes interesting 
to the automotive industry. 
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Table 3-2: Composition of the paste of the Low-Density Class-A SMC [Kec16] 


Component Trade name Quantity 
UP resin ns Premium G22-01 100 parts 
Adherent and flow aids BYK W9010 3 parts 
Styrol - 7 parts 

: Palapreg Premium G21- 
Peroxide O1LE Cure 1 part 
L&V 50%MgO Luvatol® MK35 2.77 parts 
co BE Omya Millicarb® OG 105 parts 
Filler: micro hollow glass 3M VS5500 28 parts 
spheres 


The LD Class-A SMC measurements were conducted with nearly the same 
process settings as the Class-A SMC. Due to the slightly increased thickness of 
the semi-finished product, the initial charge coverage was reduced to 36.25 % 
(290 mm x 250 mm) to gain the same final plate thickness of ht = 3 mm by 
iterative changing the initial charge dimensions. The initial charges were as- 
sembled from four single sheets which had a total thickness of ho = 9 mm. A 
representative pressure distribution for this material composition is plotted 
in Table 3-5. 
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Figure 3-4: Exemplary pressure distribution of the LD Class-A SMC for a closing speed of 


h=1mms*[Hoh17b] 


As the pressure distribution is following the same pressure pattern as that of 
the Class-A SMC, referred to in the previous Section, for the discussion of the 
pressure pattern. The slight differences of the first pressure increases of Sen- 
sor 1 to 3 (cf. Figure 3-4 and Figure 3-5) can be explained to be the slightly 
changing initial charge length and the wavy surface of the initial charge. These 
non-homogenous thickness distributions lead to a different relative gap 
height of the first pressure increase. 


3.2.3. Vinylester (VE) SMC 


The Vinylester (VE)-SMC was developed within the International Research 
Training Group (IRTG) GRK 2078 "Integrated engineering of continuous-dis- 
continuous long fiber reinforced polymer structures" by the Fraunhofer ICT. 
The objective was to develop a material formulation for semi-structural ap- 
plications. Therefore, the fiber volume content is increased, while the fillers 
are omitted. Furthermore, the material should be able to crosslink with the 
local unidirectional reinforcements during the molding process. Therefore, a 
VE resin is used. 
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In this context, the term semi-structural material is used to indicate material 
properties between the Class-A SMC and unidirectional reinforcement com- 
posites, which are used for structural applications. The semi-structural mate- 
rial needs local reinforcements or geometrical stiffeners to be used in struc- 
tural applications. [Hoh17b] 


The composition of the VE-SMC paste is given in Table 3-3. The density of the 
paste is dominated by the resin with a density of p ve-smc, paste = 1.10 g cm? as 
no fillers were used. The fiber content is set to 22.7 vol.-% (41 wt.-%) of 1 inch 
(25 mm) chopped glass fibers. Once again, the MultiStar® 272 from Johns 
Manville is used. The final semi-finished product has a density of 
p ve-smc = 1.42 g cm”. 


Table 3-3: | Composition of the paste of the semi-structural VE SMC 


Component Trade name Quantity 
VE resin Atlac XP810X 100 parts 
Adherent and flow aids BYK 9085 2 parts 
Peroxide Trigonox 117 1 part 
L&V 50%MgO Luvatol EK 100KM 4.2 parts 


During the first molding trials, the VE-SMC showed a higher flowability than 
the Class-A SMC formulations with fillers. Therefore, the initial charge cover- 
age was reduced to 20 % (160 mm x 250 mm) to investigate the viscosity over 
a longer flow length. Due to the decrease of the initial charge coverage, eight 
single SMC sheets were needed to be assembled, with a total height of 
ho = 10.1 mm. This leads to a final plate thickness of hr = 2 mm. 


The experiments were performed with the same mold temperatures, 150 °C 
for the upper and 160 °C for the lower side. The closing speed was set to 
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h=1mm s^, but also molding trials with different constant deformation rates 
of D = 0.5 st, 1.0 st and 2.0 st were performed. The deformation rate D is 
defined as the ratio between the closing speed h and the current gap 
height h: 
p = Él, (19) 
h 

To define the nominal deformation rate during the molding process, the clos- 
ing speed was prespecified over the height. During these trials, the closing 
speed varied between Amax = 20 mm s% and Amin = 1 mm s". Due to the indus- 
trial scale press and the mold size, the nominal deformation rate could not be 
held constant until the end of the filling process, because the press switches 
from velocity based control to pressure based control at the end of the com- 
pression molding process. Thus, the pressure values from sensors 5 to 7 (fast 
mold closing), or 6 to 7 (slower mold closing) could not be used for the defor- 
mation-rate-dependent experiments, since the deformation gradient de- 


creases significantly after the switching point. [Hoh17a] 


Since the experiments were performed with higher closing speeds (see Table 
3-4), the maximum press working force had to be reduced to 1 600 kN to pro- 
tect the pressure sensors. The process settings and the relevant information 
for all experiments performed with the VE-SMC are listed in Table 3-4. A rep- 
resentative pressure distribution for this SMC with a constant closing speed 
of h = 1 mm s” is plotted in Figure 3-5. 
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Table 3-4: Press settings of the VE-SMC molding trials [Hoh17a] 


: Closing Nominal de- Sensors 
Molding P y BE 
Trial Press setting speed[mm formation within set- 
s-1] rate [s] ting range 
Constant de- 
D=0.5s7 formation 2.5-5.3 0.5 1-5 
rate 
Constant de- 
D=1.0st formation 5.0- 10.5 1.0 1-4 
rate 
Constant de- 
D=2.0s7 formation 10.2-20.2 2.0 1-4 
rate 
onion Constant 
ei ' closing 1.0 0.35-0.11 1-7 
speed 
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Figure 3-5: Exemplary pressure distribution of the VE-SMC for a closing speed of 
h=1mms! 


The pressure distribution of the VE-SMC generally behaves in the same way 
as the previous SMC formulation. Nonetheless, differences in the threshold 
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pressure values and in the pressure increase in the flow phase can be ob- 
served. One clear difference is the pressure progression of Sensor 1. This has 
no SMC material specific reason, but a technical cause. It can be traced back 
to a problem with the tool. After the experiments with the Class-A SMCs, the 
short pegboard next to sensor 1 was damaged and the clearance between the 
mold sides was slightly increased. Due to this higher clearance, the tightness 
between both mold sides is not given and a slight flow occurs through this 
clearance. While closing the tool, the clearance is decreased and the fully cav- 
ity function is recovered. The loss of material can be neglected for the follow- 
ing steps, since it is very low. [Hoh17a], [Hoh17b] 


Since the pressure pattern shows the same characteristics as the Class-A 
SMCs, see Section 3.2.1 for the detailed discussion. 


3.2.4. Unsaturated Polyester Polyurethane Hybrid resin 
(UPPH) glass fiber SMC 


The UPPH glass fiber SMC was also developed within the International Re- 
search Training Group with the same aims as the semi-structural VE-SMC. As 
the VE-SMC showed a tendency to foam during the molding process, the me- 
chanical properties were not in the designated range. Therefore, the resin 
system was replaced and an unsaturated polyester polyurethane hybrid resin 
(UPPH) is used. Within the International Research Training Group, the overall 
scope is the combination of SMC with local continuous unidirectional (UD) 
reinforcements. Therefore, also UD-SMC materials were developed. Forthese 
reinforcements, the same UPPH resin is used since it has a B-stage, which en- 
ables a preforming and a rapid precuring to fixate the shape. Due to this, also 
the UPPH SMC semi-finished product can be aged in a shorter time than con- 
ventional SMCs [Büc16], [Büc17]. 


The UPPH resin based paste has a density of purrn-smc, paste = 1.08 g cm”. Its 
composition is given in Table 3-5. The fiber content of the UPPH glass fiber 
SMC is 23 vol.-% (41 wt.-%) of the same 1 inch (25 mm) chopped fibers from 
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Johns Manville (MultiStar® 272) as the previous SMC formulations. This leads 
to a density of the semi-finished product of purru-c-suc = 1.42 g cm”. 


Table 3-5: Composition of the paste of the semi-structural UPPH SMC 


Component Trade name Quantity 
UPPH resin Daron ZW 14142 100 parts 
Adherent and flow aids BYK 9085 2 parts 
Impregnation aid BYK 9076 3 parts 
Deaeration aid BYK A-530 0.5 parts 
Inhibitor pBQ 0.03 parts 
Peroxide Trigonox 117 1 part 
Isocyanate Lupranat M20R 24.2 parts 


The compression experiments were performed in the same way as the trials 
for the semi-structural VE-SMC. The initial coverage was also set to 20% 
(160 mm x 250 mm). Since the final plates should be used for mechanical test- 
ing, the final plate thickness was hr = 3 mm. Therefore, the initial charge had 
a thickness of ho = 18 mm and was assembled from 16 single sheets. 


Due to the change of the resin, the mold surface temperature had to be ad- 
justed. Therefore, the upper tool was tempered to 140 °C and the lower tool 
to 145 °C. The experiments were performed with a constant closing speed of 
h = 1 mm s" and different constant nominal deformation rates of D = 0.1 Ss“, 
0.25 st, 0.5 st and 1.0 st. Due to the increase of the initial charge volume, 
the maximum press force has been increased to 3 000 kN. During the experi- 
ments with the high deformation rates, a higher closing speed compared to 
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the previous experiments occurred due to the higher initial thickness. There- 
fore, the press pressure limit was reached too early for some experiments of 
this series. An exemplary pressure distribution of this semi-structural UPPH 


glass fiber SMC and a constant closing speed of h=1mmstis plotted in Fig- 
ure 3-6. 
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Figure 3-6: Exemplary pressure distribution of the UPPH glass fiber SMC for a closing speed 


of h=1mms? 


The pressure distribution for the UPPH glass fiber SMC shows a significantly 
different behavior compared to the other SMC formulations. At the beginning 
of the flow phase, a high pressure at Sensor 1 and Sensor 2 is recorded until 
the pressure reaches a threshold at a relative gap height of 9 mm. This thresh- 
old is at a three times higher level than the pressure at the end of the molding 
process. Compared to the UPPH carbon fiber SMC (see Section 3.2.5), this 
threshold is approx. four times higher. Since the press working pressure is 
increasing at the same time, this is not just a local effect at the locations of 
the pressure sensors, but also a global effect throughout the whole material 
at that time. After this threshold, the pressure is decreasing until it proceeds 
in the same way as the carbon SMC. From a relative gap height of approx. 3 
mm onwards, the pressure values of Sensor 1 and 2 are increasing again and 
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following the typical pressure pattern. The lower pressure values of Sensor 1 
compared to Sensor 2 can be explained once again by the damaged pegboard. 
Since the total pressure in the cavity is higher for this experiment, the mate- 
rial is still flowing and therefore, not curing occurs in the gap. Therefore, the 
pressure at Sensor 1 is lower over a larger range. Still, the amount of material 
flown through the gap is so small, that it can be neglected. [Hoh17c] 


The effect of this local pressure peak is not only appearing at the beginning of 
the flow phase, but also at the flow front along the hole molding process. 
Sensor 3 to Sensor 5 also show this behavior when the flow front is reaching 
them, but with a lower amplitude. [Hoh17c] 


3.2.5. UPPH carbon fiber SMC 


The UPPH carbon fiber SMC was developed together with the UPPH glass fiber 
SMC. Both semi-structural SMC are based on the same paste (see Table 3-5). 
In contrast to the UPPH glass fiber SMC (see 3.2.4), the fiber content was in- 
creased to 42 vol.-% (55 wt.-%) of chopped carbon fibers. The 50k carbon fi- 
ber PX 35 from Zoltek was chopped to 1 inch (25 mm) long fibers. This leads 
to a density of the semi-finished product of pupph-c-smc = 1.39 g cm”. 


The compression experiments were performed in the exact same way as the 
experiments for the UPPH glass fiber SMC. As the single carbon fiber SMC 
sheets are slightly thinner than the glass fiber sheets, the initial stack were 
assembled out of 18 single sheets. Due to the higher material costs, molding 
trials were only done for one constant closing speed of h = 1 mm stand one 
constant nominal deformation rate of D = 0.25 st. An exemplary plot of the 
pressure distribution of this SMC for a constant closing speed of h=1mms! 
is plotted in Figure 3-7. 


48 


3.3. Compressibility of SMC during flow 


mn — — Pressure Sensor 1 3:3 

150 + = — — Pressure Sensor 2 130 
sn... Pressure Sensor 3 zZ 
108 ee EEE BE | “eae Pressure Sensor 4 lose 
= — — Pressure Sensor 5 g 
E 100 Press working force | 20 Ss 
g bo 
3 £ 
3 15 = 
vo ° 
a 3 
1.0 3 
2 
a 

0.5 

T h 0.0 

0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 
Relative gap h-h; [mm] 
Figure 3-7 Emplary pressure distribution of the UPPH carbon fiber SMC for a closing speed 


of h = 1 mm s* [Hoh17b] 


The pressure distribution of this UPPH carbon fiber SMC follows the typical 
SMC pressure pattern — refer to Section 3.2.1 for the detailed discussion. The 
pressure drop of Sensor 1 after the relative gap height of 8.5 mm can be 
traced back to the problem with the damaged pegboard once again (cf. Sec- 
tion 3.2.3). 


3.3. Compressibility of SMC during flow 


In this Section, the compressibility of the different material formulations dur- 
ing the molding process is investigated. Since the SMC was assumed to be 
incompressible in most previous works (see. Section 2.3), the first material 
characterized in this Section is the VE-SMC, as the compressibility is first ob- 
served. Based on this, the other SMC material formulations are then system- 
atically examined. 
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3.3.1. Vinylester SMC formulation [Hoh17a] 


Figure 3-8 (a) and Figure 3-8 (b) show the pressure distributions of the 
VE-SMC for Sensors 1 and 4, respectively, for the three different deformation 
rates D = 0.5 s+, 1.0 st and 2.0 st. Within these plots it can be observed, that 
the relative gap heights, where the pressure is increasing, are changing with 
the different nominal deformation rates. Considering the effects described in 
the previous Section, the current averaged flow front position can be deter- 
mined and used in a correlation with the current gap height. By using the sen- 
sor positions ls, of sensor Sx, representing the position of the current melt 
front, and the current gap height h(t), the current volume of the SMC melt 
can be calculated as 


V(t) = Is, + h(t) b, (20) 


where b is the width of the mold and x the number of the sensor. Dividing the 
current volume change by the initial charge volume, the relative volume 
change ß can be calculated by 


_ 1/aV) _VEO)-Vo _ Isx h(t) lo ho 
BO) = (3) ~ VAt Ig Mg (t-te) ’ (21) 


where to represents the time when the press is touching the SMC initial 
charge for the first time. 
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Figure 3-8: Compression molding experiments of the VE-SMC: pressure over relative gap 


height for (a) Sensor 1 and (b) Sensor 4 for different nominal deformation rates 
(D = 0.5 st, 1.0 st and 2.0 st) 


The relative volume change over the real flow length is plotted in Figure 3-9 
for all constant nominal deformation rates. The first sensor reached by the 
SMC flow is Sensor 2 (S2), with a position of Isa = 165 mm, which is in the tran- 
sition between the squish and the flow phases (cf. Section 3.2.1 and 3.2.3). 
The last sensor value plotted in Figure 3-9 is the last one with a constant nom- 
inal deformation rate, before the press is switching from speed to pressure 
control. Since higher closing speeds lead to higher press forces, the switching 
point is reached earlier for higher nominal deformation rates. This can been 
recognized in the diagrams by the time, the maximum pressure threshold is 
reached at the end of the compression phase. 


The graphs plotted in Figure 3-9 show that the nominal deformation rate has 
a high effect on the volume change, where the first measuring point (Sensor 
2, | = 165 mm) does not represent the initial charge volume change, but the 
change after a short flow. After the mold is touching the material, different 
compressibility behaviors can be observed. For small nominal deformation 
rates (D = 0.5 st), a volume increase was measured, while it decreases for 
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higher nominal deformation rates. This variation can be attributed to the dif- 
ferent closing speeds and therefore different temperature and pressure dis- 
tributions during the squish phase. Higher closing speeds lead to higher ex- 
tensional viscosity (see eq. 5) and therefore to a higher pressure in the 
material. Especially during the squish (compaction) phase, this leads to a com- 
paction of the material whereby the compression of air entrapments might 
have the biggest effect. Furthermore, for small nominal deformation rates, 
the initial charge has a longer contact with the lower molding tool (1 s, +10% 
contact time), so the temperature is slightly higher. This increases the release 
of styrene and therefore a volumetric increase can occur. Additionally the 
void and gas entrapments content might change, like shown by Comte et al. 
[Com06] and Lee et al. [Le08]. 


At the beginning of the flow phase, the tendency from the squish (compac- 
tion) phase is the same, but approximately after a flow length of 240 mm (at 
Sensor 4, | = 392 mm), the volume change reaches a threshold at £ = 0 s. 
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Figure 3-9: Comparison of the relative volume change £ of the VE-SMC between measure- 


ment and modelling for three nominal deformation rates D = 0.5 st, 1.0 stand 
2.0 s? during the molding process [Hoh17a] 
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This behavior clearly shows that this semi-structural SMC material cannot be 
assumed as incompressible. The most common approach to describe the re- 
sulting volume change induced by a pressure is to use the bulk modulus. Since 
the bulk viscosity can only describe a volumetric decrease due to a pressure 
increase, this method is not suitable for this case. Another approach to de- 
scribe the volume change is to use an empirical model to predict the changes 
during the flow. By evaluating the experimental data (compare Figure 3-8), it 
could be observed that the relative volume change p depends on the defor- 
mation rate and the flow length. Combining this with a power-law estimation 
and an empirical evaluation, the following empirical equation is evolved by 
fitting different equation terms: 


max” "incom (t) k eT = 
b(t, DC) = a, (Te) (D, = DW) - on” , (22) 
ldef 


with laef = 1m, the incompressible flow front position lincomp defined by the 
mass balance equation as 


loh 
Lincomp(t) = ae (23) 


with the power-law parameters k and |, the fitted volume change factor a, 
and a, and the deformation rate D., where the compressibility is changing 
from a volume decrease to a volume increase. The comparison between the 
measurements and the empirical model are plotted in Figure 3-9. 


3.3.2. Class-A SMC formulations 


After showing the compressibility of the VE-SMC, also the other SMC formu- 
lations are evaluated. In Figure 3-10, the averaged experimental relative vol- 
ume changes ß are plotted for both Class-A materials and a constant closing 
speed of h = 1 mm s^. As the relative volume change is nearly zero over the 
whole flow length, the Class-A SMC formulations can be assumed to be in- 
compressible. 
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Figure 3-10: Averaged relative volume change £ from the Class-A SMCs experiments for a 


constant closing speed of h=1mms? 


3.3.3. UPPH SMC formulations 


For the UPPH glass fiber SMC, the evaluation of the relative volume change 
was only possible for the constant deformation rate of D = 0.25 st. For the 
other test series, the fingering effect was too high. The fingering effect de- 
scribes the statistically forming of faster flowing “fingers”, due to local varia- 
tion of the viscosity [Yad11]. This effect can be recognized in Figure 3-6 for 
the sensors 3, 4 and 5: While the flow front is reaching the sensor, a first small 
pressure increase occurs. When the finger slows down, the pressure stays 
constant, until the average flow front reaches the sensor position. Then the 
pressure is increasing again. This effect can also be observed for the VE-SMC 
(cf. Figure 3-5) at Sensor 4, when a pressure is recorded for the first time at 
2.5 s. The threshold of approx. 3.0 bar is the pressure of the small finger, while 
the normal pressure increase due to the rheological pressure starts at approx. 
2.2 s. For the UPPH glass fiber SMC, this effect is much higher and overlaid by 
other effects (cf. Section 3.7). Therefore, the determination of the exact time 
at which the flow front reaches the sensor has a high inaccuracy. Since only 
the experimental series for the deformation rate of D = 0.25 s+ has enough 
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repetitions to average the inaccuracies, the relative volume change is only 
plotted for this setting (see Figure 3-11). In the same plot, the relative volume 
change of the UPPH carbon fiber SMC is given. Based on this experimental 
data, the UPPH SMC formulations have to be considered as compressible. 
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Figure 3-11: Relative volume change £ from the experiments of the glass and carbon UPPH 


SMC for a nominal deformation rate of D = 0.25 st 


3.4. Compressible Rheological 2D Shell-Model 
[Hoh17a] 


After modelling the compressibility, a compressible rheological shell model is 
developed to take into account the new insights from Section 3.3. The basic 
rheology model used in this Section was developed by Dumont et al. [Dum03], 
[Dum07], [Org12] and will be modified to fulfill the compressibility condition. 
For incompressible SMC formulations, the same model can be used by using 
a constant relative volume change ß = 0. 


From previous studies, Dumont et al [Dum07] made the following assump- 
tions for a plug flow model: 
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SMC is a one-phase continuum medium. This means, that the paste (resin 
and filler) and the fibers have the same velocity and no fiber-matrix sep- 
aration will occur 

The temperature in the material is practically isothermal [Bar79], except 
in the thin resin layer near the mold 

There will be no curing during the compression molding. As shown previ- 
ously, the filling time is less than 15 s and together with the isothermal 
assumption, the reaction does not start during mold filling [Abr00] 

No external or acceleration forces will be considered, as they are very 
small relative to the viscous forces due to the deformation. Therefore the 
stress tensor is defined as 


o = -pö+oy(D), (24) 


where p is the hydrostatic pressure, oy the viscous stress tensor and 
D the deformation rate tensor. Therefore, the first momentum bal- 
ance equation is defined as 


div(o) = 0 (25) 


Based on Section 3.3 the volumetric change is assumed spatially constant 
and isotropic [Gui12] 


By using the coordinate system introduced in Figure 3-1 and using the middle 


plane of the current plate height as x3 = O and the x3 boundary conditions at 


the mold walls, Dumont et al. [Dum07] derived the following formulation for 


the strain rate : 


1 
v1 a (v12 + v21) 0 
1 
3 (v12 + V21) V2,2 0], (26) 
0 0 2 
h 


where vij is defined as: 
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(27) 


As shown in the previous Section, the material is considered as compressible. 
This leads to the following mass balance equation: 


div(@) = -2= Elt) = B(t), (28) 


with the velocity field ©, the density p, its change p and the volumetric strain 
rate &,. 


Due to the plug-flow behavior, the 3D problem can be reduced to a 2D shell 
description in the e1-e2-plane. Due to the fact that the shear components D,3 
(y E {1,2}) of the strain rate tensor are constrained to zero, arbitrary reaction 
terms (T,3) must be added [Dum07]. Therefore, the stress tensor ø from eq. 
24 can be rewritten as: 


o= -pöH+T,s(e, 8 e +83 @ ey) +oay, ye {1,2}. (29) 


To gain the momentum balance equation (eq. (25)) for the shell-approach, 
the equation is integrated over the thickness h of the SMC [Dum07]. There- 
fore, it can be written as: 


-gradP + (T+ + r-) + div’ = 6, (30) 


~My 


where the symbol “” indicates the 2D unknown field in e1-e2-direction of the 


3D tensor, with the thickness averaged pressure P and the thickness averaged 
in-plane viscous stress È”, defined as: 


pati? pax (31) 
h My, p 39 


X= if v2 BY dx (32) 
= h My, 30 
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and the tangential stress vectors T+ and T- between the SMC and the upper 
(x3 = h/2) and the lower (x3 = -h/2) part of the mold, respectively: 


= N LL, h 
(+7). a= [Ts] 2, y € {1,2}. (33) 


By using the assumption that no external or acceleration forces will be con- 
sidered (eq. (25)), and that the deformation rate is constant in the material at 
one point of time, which leads to 

aa’ _ 


—_=0, (34) 


dx; 
the thickness averaged mass balance equation (eq. (30)) results in: 


dP ıd(T*+T7) =Ü 
dxı h dxı = , (35) 


for the rheological tool flow. By using the power-law viscosity model pro- 
posed by Dumont et al. [Dum03] and the power-law friction stress model of 
Abrams and Castro [Abr03] modified to the 2D approach by Dumont et al. 
[Dum07], the averaged pressure P(xı) can be solved by integrating these 
equations along with the boundary condition at the flow front. The boundary 
condition for this integration differs from the incompressible one, since the 
compressible flow front (x1) is defined by equation (21), where Is = x1: 


xl) = 22 (ate BO +1). (36) 


This leads to the following equation for the thickness-averaged pressure with 
an initial charge sample at the left side: 
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hl 


I ae 


Pa) = -mla i 


1+H/ D 


( 2a ) h 
m+1/ vin-t p2 


where „sis the corresponding extensional viscosity (compare Barnes 


A 
h 


| a - (2% (at BC) + D] , (37) 


h(t) 


[Bar00]) under a plane strain rate of Dy = 1.0 st, H is a rheological function 
(compare Dumont el al. [Dum07]), n and m are the power-law exponents for 
the viscosity of the SMC and for the rheology of the boundary layer, respec- 
tively, and A is the hydrodynamic friction coefficient. 


Using the global form of the momentum balance equation, averaging it over 
the volume and inserting eq. (37) (compare Dumont et al. [Dum07]), leads to 
the following equation for the volumetric averaged axial stress (033): 

a er 
nh 


h 


h 


h 


1 


(033) = ps pn-T 


2A \ A AAN floh m+1 
beim) (ion) > (38) 
By using eq. (30), the local thickness averaged axial stress X33 at the mold 
surface (x3 = +/- h/2) can be determined by 


1 


n-i 
Do 


h 


233 (%1) = Nps h 


HT en-eo)". e 


mti vo 


These equations can be seen as the sum ofthe bulk rheology stress (first term) 
and the friction stress (second term). 
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3.5. Material Characterization with the 
Compressible Rheological Model 


To determine the rheological parameters of the core region and the friction 
parameters of the lubrication layer, a two-step approach is used. 


First, the hydrodynamic friction coefficient A and the corresponding power 
law coefficient m are identified using the pressure difference between two 
sensors at a constant closing speed. By using equation (39) at the mold sur- 
face (x3 = +h/2), this pressure difference can be expressed as: 


033(Xır) — 033(X;) = — am i = kurt - xp]. (40) 
m+1 \vo h 

In the second step, the viscosity power-law coefficient n is fitted to the exper- 
imental pressure results at different nominal deformation rates using equa- 
tion (39). To gain a representable strain-rate range, the results for the con- 
stant nominal deformation rates and the constant closing speed have to be 
used [Hoh17a]. Based on this coefficient, the reference extensional viscosity 
Nps can be determined. 


In the following Section, this procedure is applied to the test series of the dif- 
ferent SMC formulations. Furthermore, the resulting material parameters are 
presented. The comparison of all material parameters are given and discussed 
in Section 3.6. 


3.5.1. Class-A SMC formulations 


For the Class-A SMC formulations, molding trials only with a constant closing 
speed were performed. Since the hydrodynamic friction coefficient charac- 
terization dependents on the closing speed (cf. Section 3.5) and not on the 
deformation rate, the variation over the process range is relatively low. The 
experiments with a constant deformation rate would have led to a higher var- 
iation within the molding series. Based on the available data, a friction power- 
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law coefficient of m = 0.6 is determined for both materials. The hydrodynamic 
friction coefficient exhibits a small variation between the formulations. The 
normal Class-A SMC has a hydrodynamic friction coefficient of 
Aciass-A = 1.8 MNs m”, while the LD Class-A SMC has a slightly higher hydrody- 
namic friction coefficient of An class-a = 2.0 MNs m”. 


For the LD Class-A SMC, additionally experiments were performed with paste 
without fibers. By using the same method, a friction power-law coefficient of 
m=0.6 and a reduced hydrodynamic friction coefficient of 
ALp Class-A, paste = 0.85 MNs m” are determined. 


After the determination of the friction parameters of the lubrication layer, 
the reference extensional viscosity is calculated. In a first step, the rheological 
stress ov33 is plotted over the deformation rate D to determine the rheological 
power-law exponent n. Therefore, the reference deformation rate of 
Do= 1 stis used. This correlation is shown in Figure 3-12. For the Class-A SMCs 
the exponents are identified with nciass-asmcs = 0.44. For the LD Class-A SMC, a 
coefficient of nıo ciass-A smc, paste = 0.58 is evaluated. These two coefficients for 
the SMC and the paste are the same results as determined by Dumont et al. 
[Dum03]. [Hoh17b] 


Using these rheological power-law coefficients, the reference viscosity for the 
Class-A formulations can be determined as Nps class-asmcs = 4.10 MPa s and for 
the paste as Nps LD Class-A SMCs, paste = 0.75 MPa s. 
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Figure 3-12: Fitting of the viscosity power-law exponent n with the results for the constant 
closing speed (A = 1 mm s?) of the Class-A, LD Class-A SMC and the LD Class-A 
SMC paste 


3.5.2. Vinylester SMC formulation 


For the VE-SMC characterization, a high number of experiments are available, 
covering a large range of process settings. Therefore, the material data could 
be determined over a large closing speed and deformation rate range with 
many data points (cf. Section 3.2.3). 


By using the pressure difference expressed by eq. (40), a friction power-law 
coefficient of m = 0.6 and an averaged hydrodynamic friction coefficient of 
Ave-smc = 2.3 MNs m”? is determined. These values are close to those presented 
by Dumont et al. [Dum07]. 


After determining the friction coefficients, the viscosity and the viscosity 
power-law coefficient can be computed. These fitted results, together with 
the trend line, are plotted in Figure 3-13. By using the reference deformation 
rate of Do=1.0s7 in eq. (38), the best viscosity power-law exponent is 
Nve-smc = 0.27, with a small variation over the different experiments. This 
value corresponds to the ones presented by Lee et al. [Lee81]. [Hoh17a] 
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Figure 3-13: Fitting of the viscosity power-law exponent n with the results for the constant 


nominal deformation rates (D = 0.5 st, 1.0 s? and 2.0 st) and the constant clos- 
ing speed (A = 1 mm s?) of the VE-SMC formulation [Hoh17a] 


3.5.3. UPPH SMC formulations [Hoh17b] 


Due to the pressure distribution with the pressure peaks at the flow front of 
the UPPH glass fiber SMC, compared to the UPPH carbon fiber SMC (cf. Sec- 
tion 3.2.4 and 3.2.5), the characterization of the rheological parameters is 
only possible for the carbon fiber SMC. The presented compressible rheolog- 
ical model (cf. 3.4) cannot describe material property changes during the 
flow, which must be the reason for the pressure peaks at the flow front. To 
give a first description of these material property changes, a general model- 
ling approach for the glass fiber SMC is presented in Section 3.7. 


During the evaluation of the hydrodynamic friction coefficient for the UPPH 
carbon fiber SMC, it becomes apparent, that it is not possible to determine a 
constant value for the velocity controlled nor for the deformation-rate-con- 
trolled experiments. Comparing the three diagrams for three different time 
steps during the molding process in Figure 3-14, it can be recognized, that the 
coefficient has two thresholds. One is in the area of the initial charge (flow 
front smaller than 0.16 m) and the other one is in the flow area. This becomes 
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even clearer when taking the small leakage and therefore the pressure de- 
crease of Sensor 1 into account. When the flow front reaches Sensors 4 and 
5, the pressure at Sensor 1 is underestimated and therefore also the pressure 
difference. This causes the decrease of the hydrodynamic friction coefficient 
Aatx=0.035 m. 
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Figure 3-14: 
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Progression of the hydrodynamic friction coefficient A over time and location 
for the UPPH carbon fiber SMC and a constant closing speed of h = 1 mm st 
[Hoh17b] 
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The reason for these two different hydrodynamic friction coefficients cannot 
be related to the temperature nor to the thickness of the lubrication layer, as 
this would lead to a smaller value in the area of the initial charge. Since this 
effect does not appear for the glass fiber reinforced materials, it could be 


caused by the interaction of the carbon fibers with themselves or with the 
resin. 


By plotting the axial rheological stress over the deformation rate, a best vis- 


cosity power-law exponent of nureH c- smc = 0.27 could be fitted (see Figure 
3-15). 
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Figure 3-15: Fitting of the viscosity power-law exponent n with the results for the constant 


nominal deformation rates (D = 0.25 s?) and the constant closing speed (Å = 1 
mm s!) of the UPPH C-SMC formulation 


As detected for the friction characterization, the UPPH carbon fiber SMC 
shows a variation also for the extensional viscosity in the initial charge and 
the flow area. The averaged value for the extensional viscosity in the IC area 
is determined to npsic = 3.65 MPa s, while it is nps flow = 2.75 MPa s in the flow 
area. Again, this could not be ascribed to temperature distribution, lubrica- 
tion layer thickness or other process parameters. Considering both variations 
for the friction and viscosity stresses, the anisotropy of the viscosity due to 
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fiber orientations and fiber volume changes could be an explanation. Both 
effects were numerically investigated by Tsigkopoulos et al. [Tsi09] and Ber- 
töti et al. [Ber17]. They developed numerical anisotropic viscosity ap- 
proaches, which describe the viscosity as a function of the fiber orientation. 
This correlates to the experimental observation that the difference between 
the initial charge area and the flow area depends on the closing speed (see 
Figure 3-16). The change ofthe closing speed leads to a different deformation 
gradient and therefore to different fiber orientation distributions. The reason 
why this effect only appears for the carbon fiber reinforced SMC might be the 
fiber bundle structure. For glass fibers, the bundles stay in a more or less 
round shape, while for carbon fibers, the bundles spread open to a flake-like 
structure (see Figure 3-17). This effect is intensified by the different bundle 
sizes, since the UPPH carbon fiber SMC is manufactured with 50k rovings, 
while 12k rovings are used for the UPPH glass fiber SMC. Therefore, other 
fiber interactions and fiber reorientation mechanisms take place. Another 
possible explanation for this behavior could be based on the different thick- 
ening mechanisms of Class-A respectively VE SMC compared to the UPPH 
SMC. The first two are thickened by an alkaline earth metal oxide (magnesium 
oxide). This thickening mechanism is investigated since 1970 [Van70] but still 
not understood satisfyingly [Kec16]. The UPPH SMC is thickened with isocya- 
nate. Hereby, the hydroxyl and carboxyl groups of the UP form covalent 
bonds with the isocyanate and consequently increase the extensional viscos- 
ity of the paste [Sai96]. As the used isocyanate shows a reactivity of 2.7, a 
weak three-dimensional network will be the result of the cross-linking during 
maturation. Compression molding of the SMC should cause a partially break 
of this network at a certain strain. This would consequently lower the exten- 
sional viscosity of the SMC during molding. 
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Figure 3-16: Relative change of the extensional viscosity and hydrodynamic friction between 


the initial charge area and the flow area of the UPPH carbon fiber SMC 
[Hoh17b] 


Figure 3-17: Comparison of the fiber bundle structure in the glass fiber reinforced VE SMC 
(left) and the UPPH carbon fiber SMC (right) 
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3.6. Effects ofthe material type on the rheological 
behavior and model parameters 


Based on the variations in the SMC material composition, effects like com- 
pressibility and variations in the initial charge area and the flow area are ob- 
served. The Class-A SMC with fillers and a relatively low glass fiber volume 
content showed the typical rheological characteristic with similar material pa- 
rameters to those presented by Dumont et al. [Dum03]. By increasing the fi- 
ber volume content and omitting the fillers, a substantial compressibility 
could be observed. The change from glass to carbon fibers leads to material 
behavior, which shows the same compressibility, but a significant change be- 
tween the initial charge area and the flow area. This could be related to the 
structure and shape of the fiber bundles or of the matrix or both. While show- 
ing their specific differences in rheological behavior, all SMC formulations, ex- 
cept for the UPPH glass fiber SMC, can be modeled with the same basic model 
presented by Dumont et al [Dum03], [Gui12], [Fer17]. Within this model, only 
small modifications regarding compressibility and local dependencies (IC and 
flow region) of the viscosity and friction parameters due to local variation of 
fiber orientations according Section 3.4 and 3.5 are necessary. 


Taking these information and characterization results into account, it is pos- 
sible to choose a suitable modelling and characterization approach for a spe- 
cific SMC. This helps to simplify the characterization procedure and to set the 
results into the correct context. 


Furthermore, the power law coefficient for the friction term m can be as- 
sumed to be constant, independent of the material composition. The exten- 
sional viscosity coefficient n can be considered as constant depending on the 
SMC type. For SMC formulations with a high filler content it is 0.44, while the 
pure resin mixture of such an SMC has a coefficient of 0.58. By omitting the 
fillers for semi-structural SMC formulations, a coefficient of 0.27 can be used. 
Based on this, only the extensional viscosity n and the hydrodynamic friction 
coefficient needs to be determined for new SMC formulations. 
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Table 3-6: Comparison of rheological material data for a closing speed of h=1mm s! (comp. 


[Hoh17b]) 
N n A Hydrodynamic fric- 
Extensional viscosity a ; 
aai tion coefficientA ath m 
nesatD=15 u 
=1mms 
Class-A 4.10 MPa s 0.44 | 1.80 MN s m° 0.6 
LD Class-A | 4.10 MPa s 0.44 | 2.00 MN s m” 0.6 
ED Class ep 0.58 | 0.85 MN sm? 0.6 
resin 
VE SMC 12.00 MPa s 0.27 | 0.80 MN sm” 0.6 
3.65 MPa s 1.23 MN s m° 
(IC area) (IC area) 
OPRSNICE 2.75 MPa s 07 0.60 MN sm” ne 
(flow area) (flow area) 


3.7. Alternative flow kinetic of the UPPH glass fiber 
SMC [Hoh17c] 


In Section 3.2.4, 3.3.3 and 3.5.3 the different behavior of the UPPH carbon 
fiber SMC, compared to the other SMC formulations, are discussed. To iden- 
tify the reasons for this behavior, an experimental investigation based on op- 
tical and micro tomography scans were done. Based on these results, an al- 
ternative flow description is presented. 


3.7.1. Experimental investigation 
Optical and pressure based analyses 
To identify the reason for the pressure peaks at the flow front of the UPPH 


glass fiber SMC, a series of short shots was conducted with two repetitions 
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for every step. During the first evaluation of these plates, three different ar- 
eas are observed by optical assessment. At the flow front and at the opposite 
end of the plate (where the initial charge was placed), areas with a random- 
ized oriented fiber distribution and complex curvature of the fibers at the sur- 
face are observed (Figure 3-18). In the center area between the other two 
areas, a quasi-isotropic fiber orientation with mainly straight fibers is ob- 
served at the surface. 


CT Sample 2 
(center section) 
= 


Flow 
direction 


Figure 3-18: Gray-scale picture with high contrast of a UPPH glass fiber SMC short shot and 


CT Sample positions 


To identify which of three regions is responsible for the pressure peaks, the 
local lengths of the peak are evaluated. Therefore, the point with the first 
significant pressure decreases after the pressure peak is used. By using the 
press data with the current gap height, the position of the current flow front 
is calculated [Hoh17a]. Using this flow length information and the position of 
the sensor with the pressure peak, the length of the primary responsible area 
could be calculated. Furthermore, the plates were gauged by using visual cri- 
teria according to Figure 3-18. Therefore, the three areas were measured. The 
results from the pressure evaluation and the length of the swirled flow front 
region (vortex area) are plotted in Figure 3-19 for both plates (repetitions) of 
the short-shots study with 4.5mm thickness. The optical method leads to only 
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one measuring point, but with the pressure evaluation method, a tracking of 
the flow front region could be done. As the correlation between the optical 
and the pressure evaluation seem to correspond, the flow front region with 
its highly swirled fiber orientation at the surface seems to be responsible for 
the pressure peak. 
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Figure 3-19: Length of the vortex area at the flow front: optical measured and extracted 
from pressure distribution 


Micro tomography scans 


To understand what is happening during the flow phase and why the three 
regions develop, a series of micro tomography scans were conducted by the 
Institute for Applied Materials (IAM-WK) at the Karlsruhe Institute of Tech- 
nology (KIT). Therefore, four samples of 20 mm x 20 mm, two from the flow 
front region and one each from the center and initial charge region were ex- 
cised from each 4.5 mm thick plate. The positions of these samples are indi- 
cated in Figure 3-18. 
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Comparing the scans from the flow front (Figure 3-20) and the middle region 
(Figure 3-21), the visual differences between these two regions can be con- 
firmed. The flow front sample shows fibers, which are highly swirled. Further- 
more, the original fiber bundles have been separated into smaller bundles or 
even single fibers at the flow front. In contrast to these vortex regions, the 
center region shows the typical fiber orientation originating from the manu- 
facturing process. Especially in the cross-section, the randomly orientated fi- 
ber bundles can be recognized. At the surface, some separation and reorien- 
tation effect of the fiber bundles appeared. 


Figure 3-20: Micro tomography scan with color-coding of the fiber orientation direction of 
the flow front (sample 1): view on top of the sample (left) and sliced in the mid- 
dle (right). 
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Figure 3-21: Micro tomography scan with color-coding of the fiber orientation direction of 
the sample from the center region (sample 2): view on top of the sample (left) 
and sliced in the middle (right). 
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3.7.2. Model Discussion 


The typical flow model, which is also used for the rheological characterization 
of the carbon fiber SMC [Hoh17a], divides the material in two regions in thick- 
ness direction. Since the SMC initial charge is placed room-tempered into the 
hot mold, the outer layer of the material is heating up rapidly. Owing to the 
relatively low thermal conductivity, this layer has a varying thickness, in the 
area of only a few micrometer like shown by Shokrieh and Mosalmani 
[Sho14]. This layer contains no fibers and, because of the high temperature, 
the viscosity decreases dramatically and the layer functions as a lubrication 
layer. The charge core region stays relatively cold during the molding process. 
Due to these two different regions, the lubrication layer is exposed to a high 
shear rate, while the core charge is loaded by normal forces. Since SMC is a 
compressible, like shown in Section 3.3, and a Non-Newtonian fluid [Kot98], 
[Gib99], the shear and the extensional viscosity are not coupled by the Trou- 
ton ratio [Tro1906] and have to be determined separately. 


The extensional-based deformation of the core region leads to a material 
flow, where the core charge stays in the middle of the plate thickness and no 
mixing in the thickness direction occurs. Therefore, the measured pressure 
can be divided into a rheological stress (from the core charge region) and the 
hydrodynamic friction stress (from the lubrication layer) (cf. eq. (1)). 
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© Lubrication Layer 


Core region 


Flow 
direction 


Figure 3-22: Schematically cross-Section of SMC material during the molding process 


In contrast to this typical behavior of the carbon fiber SMC, the glass fiber 
SMC must have a different flow behavior to explain the pressure distribution. 
To identify this flow behavior, the evolutions of the three regions (vortex re- 
gion adjacent to the mold walls and aligned region in the core region) over 
the molding process is analyzed using the short-shot study. Taking into ac- 
count the micro tomography scans it seems that the core charge is moving in 
flow direction. Therefore, the lubrication layer or more precise in this case, 
the complete sheet at the mold surface is flowing into the vortex regions (see 
Figure 3-23). Furthermore, the core charge region is mainly transported in 
flow direction, while showing low elongation. In such a system, the main flow 
of the material would be at the mold surfaces, where the higher temperature 
leads to a lower viscosity. 


The reason for the pressure peak could not be identified yet. It can be ex- 
cluded, that the material is speeding up so much, that the hydrodynamic pres- 
sure becomes relevant. With a flow front speed of v = 0.1 ms“ after 4 s of the 
molding process and a 10 times higher speed in the lubrication layer, this 
would only lead to an increase of 1.0 bar. It is more realistic, that the fiber- 
fiber-interaction in the vortex regions leads to an increase of the homoge- 
nized shear viscosity and therefore to higher stresses. Furthermore, the fiber 
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bundle separation and the curved fibers could cause an anisotropic viscosity 
shown by Bertöti and Böhlke [Ber17] and Tsigkopoulos et al. [Tsi09], which is 
higher than the isotropic viscosity. New methods to investigate this approach 
will be designed and conducted in the future within the GRK2078. 
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Figure 3-23: Schematically cross-Section of the probable flow of glass fiber UPPH SMC dur- 
ing the molding process 


These two different flow behaviors are caused by the fiber type, as the paste 
is the same for both UPPH SMC materials. Therefore, the thermal conductivity 
and consequently the temperature distribution seem to be relevant. As the 
UPPH is still a new resin under research, the material parameters, and the 
material property effects during the compression molding process have to be 
investigated further. Furthermore, the ratio of the fiber stiffness to the matrix 
viscosity might have a significant effect [Pér15]. However, no reason for these 
two different flow regimes (see Figure 3-22 and Figure 3-23) of the glass and 
carbon fiber UPPH SMC could be verified so far. The UPPH resin will be further 
investigated within the International Research Training Group “Integrated en- 
gineering of continuous-discontinuous long fiber-reinforced polymer struc- 
tures” and new results will be available within the next years. 
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3.8. Summary on rheological measurements 


Chapter 3 “Rheological measurements” shows the difference of the flow be- 
havior of different SMC formulations. For this purpose, a rheological in-mold 
tool was developed. This rheological tool has seven pressure sensors aligned 
in flow direction. By recording these pressure data and the press data simul- 
taneously, a characterization of the SMC formulations is possible. 


Standard Class-A material formulations show a flow regime and a resulting 
pressure distribution, which could be described by the models developed by 
Dumont et al. [Dum03], [Dum07], [Gui12]. Using this model, the two Class-A 
formulations presented in Section 3.2 are characterized. 


To capture new application areas for SMC as a semi-structural material, the 
fiber content of SMC is increased while omitting the fillers. By changing the 
composition in this manner, it was shown, that these semi-structural SMCs 
could not be assumed to be incompressible anymore. Therefore, compressi- 
bility was added to the previous model by adding the relative volume change 
as a specific material property. This model was validated by using the semi- 
structural VE SMC. 


To combine the semi-structural SMC with local unidirectional reinforcements, 
a new unsaturated polyester polyurethane hybrid resin (UPPH) was used for 
both components, the continuous and the discontinuous fiber reinforced ma- 
terials. This UPPH resin has a B-stage to fix the preformed shape of the local 
reinforcements. The B-staging is also effecting the flow behavior of the dis- 
continuous fiber reinforced SMC material. It was shown, that the carbon fiber 
SMC has different material properties in the initial charge area and the flow 
area, but could be described by the compressible rheological model. By 
changing the fiber type from carbon to glass fiber, this UPPH SMC is changing 
its flow behavior. It was presented, that the flow front has high-pressure 
peaks, while the remaining areas follow the normal typical pressure pattern. 
Using these pressure distributions and additional CT scans, a new flow pattern 
was proposed. This characterization can be used as a basic assumption to 


77 


3. Rheological measurements and rheological modelling 


build up an alternative stress model to Dumont et al., which can describe SMC 
formulations, which do not have a distinctive bulk flow. 
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In the first section of this chapter, the current state of the art process simula- 
tion is evaluated by using the pressure distributions of the rheological meas- 
urements of Chapter 3. Due to the inaccuracy of the obtained pressure pre- 
dictions of the state of art software, the Coupled-Eulerian-Lagrangian (CEL) 
approach is used to implement a new SMC-specific material model. There- 
fore, section 4.2 is focusing on the basic modelling and the implementation 
of these models. Furthermore, a first consideration of homogeneous local re- 
inforcements in the compression molding process is given. In section 4.3, the 
validation with the pressure results from section 3.2 is performed, as well as 
a quantitative validation with the IRTG reference structure. 


4.1. Application and Evaluation of state of art 
commercial compression molding software 


To evaluate the state of art commercial compression molding software tools, 
the characterization pressure distribution of the Class-A SMC (cf. Section 
3.2.1) is used. Therefore, the mold is modeled with the exact sensor position 
within Autodesk® Moldflow® to compare the pressure results at the exact 
same position. Due to the equation of state formulation in Moldflow, the SMC 
material is compressed at the beginning of the molding process. To compen- 
sate this compressibility, the initial charge height is increased slightly from 
8.4 mm to 9.0 mm. The model setup is shown in Figure 4-1. 
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0302 


Scale (800 mm) 


Figure 4-1: Compression molding model with Moldflow, green: final cavity; red: initial 
charge position; black points: sensor positions 


Since Moldflow has its historical origin in the injection molding simulation of 
thermoplastic materials, the shear- and extensional viscosity are coupled by 
the Trouton ratio [Tro1906] for incompressible Newtonian fluids by 


Nshear = 1 . (41) 


Next 3 


Therefore, only the shear viscosity is calculated by using the Castro-Macosko 
[Cas82] reactive viscosity model 


(Cy +C2@) 
PER no(T) ag ; 
n(a,T,y) = (mon FT (m on (5) , with (42) 
= Tp 
no(T) = B exp (2), (43) 


which is an extended Cross Williams-Landel-Ferry (Cross-WLF) model for re- 
active resins. In these equations, t*, B, Tą, C4, C2 and a, are data-fitted coef- 
ficients, æ is the current degree of cure in the range between 0 and 1. To fit 
the material parameters of the reactive viscosity model, the pressure distri- 
butions from the Class-A SMC characterization (see Section 3.2.1) are used. 
Using the best fit over all seven pressure sensors, the pressure distribution 
shown in Figure 4-2 is predicted by the commercial tool. At the beginning of 
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the compression molding process, the simulated curves (grey) do not have 
the characteristic pressure increase, which is followed by the thresholds. Dur- 
ing the flow phase, the simulated pressure shows the same pattern as the 
experimental data, but it is overestimating the pressure towards the end of 
the compression mold process. Due to this overestimation, the maximum 
press force is reached too early and therefore, the pressure decrease due to 
the reduced closing speed begins at a higher relative gap height. 
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Figure 4-2: Comparison between the pressure distribution from the experiments (black) 
and the commercial simulation tool based on the fitted reactive viscosity model 
(grey) 


Another important criterion for the comparison is the fiber orientation. In Fig- 
ure 4-3, the degree of alignment (first eigenvalue) of the fiber orientation ten- 
sor A is plotted for the final plate thickness on the top and bottom surface as 
well as in a center cutting plane. The fiber orientation of the initial charge is 
defined as quasi-isotropic, which is a typical assumption for SMC. Therefore, 
the first and second eigenvalues (equal A11 and A22) are set to 0.5. To predict 
the fiber re-orientation during the molding process, the Folgar-Tucker model 
with the ARD-RSC extension (cf. Section 2.4) is used. In the fiber orientation 
plots, which are distributed over the thickness via three Sections, two effects 
can be observed. The first one is that the fiber orientation on the two surfaces 
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of the plate is not changing over the flow length. Only near the left tool wall, 
where the initial charge was placed, a slight decrease of the alignment can be 
observed. Despite the initially identical fiber orientation, both surfaces show 
a slightly different fiber orientation. The second effect is that the fiber orien- 
tation in the center plane shows a different behavior. Near the short tool side, 
where the initial charge was placed, the maximal fiber alignment can be ob- 
served. In the center, the gradient is opposed to the one at the surfaces. The 
fiber alignment decreases over the flow length. 


ap! cT2 
Oo Oo 
Figure 4-3: Alignment of the fiber orientations predicted by the commercial simulation tool 


by starting with a quasi-isotropic fiber orientation distribution in the initial 
charge: top: Top surface view middle: Cutting surface in the middle plane bot- 
tom: bottom surface view 
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To validate the fiber orientation distribution, two CT-scans were made by the 
Institute for Applied Materials (IAM-WK) [Bry14] to evaluate the distribution 
through the thickness. Moreover, a series of tensile tests in and perpendicular 
to the flow direction has been performed in previous work and can be used 
to evaluate the anisotropy over the flow length. 


The two CT-scans were extracted at approximately half the plate’s length. The 
exact position is given in Figure 4-3. To compare the fiber orientation with the 
predicted ones, the three fiber orientation tensor components A11, A22, A12 
are compared, since they describe the 2D fiber orientation in the plate plane 
(see Figure 2-9). These distributions over the thickness are given in Figure 4-4. 
As the sample size is smaller than the fiber length, a small deviation between 
the fiber orientations can be accepted, since SMC shows local statistically dis- 
tributed scatterings. Therefore, the prediction of A22 and A12 fiber orienta- 
tion distribution components seems to match, while A11 is predicted consid- 
erably smaller than measured by CT. Taking into account the form of the fiber 
orientation tensor and the plate’s thickness to fiber length ratio, the sum of 
A11 and A22 should be equal to one. Since no fibers could be orientated in 
thickness direction, the A33 component has to be zero. However, it can be 
observed, that this criterion is not fulfilled for the commercial process simu- 
lation tool. It has constant A33 component of approximately 0.2, which is 
could not be detected in reality. 
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Figure 4-4: Comparison between the fiber orientation distribution from the simulation and 


CT-scans (from [Bry14]): CT-Position 1 at x = 339 mm, y = 200 mm (top) and CT- 
Position 2 at x = 431 mm, y = 200 mm (bottom) 


To evaluate the anisotropy over the flow length, tensile tests are used, as they 
are less elaborate than CT measurements. The resulting Young’s moduli in 
both in-plane directions are given in Figure 4-6. This figure clearly shows a 
uniform anisotropy over the flow length, with some small deviations and 
slight effects close to the wall at the initial charge position. In contrast, the 
fiber orientation distribution from the commercial process simulation tool 
shows an anisotropy change during the molding process along the flow length 
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(cf. Figure 4-4). This anisotropy change is relatively high where the initial 
charge was placed, but also present in the flow area, where a smaller fiber 
reorientation occurs. These orientation effects can be explained by the shear 
flow domination in the commercial tool. Figure 4-5 shows two cross-sections 
of the flow front during the compression molding process. The shear rate and 
the velocity in flow direction are showing the typical thermoplastic fountain 
flow behavior, which has nothing in common with the SMC typical plug-flow. 


00900 


Figure 4-5: Shear rate (top) and velocity in flow direction (bottom) at the flow front during 
the compression molding simulation with the commercial simulation tool 


Considering the fiber reorientation of the combined ARD-RSC Folgar-Tucker 
model (Section 2.4, eq. (11)), it can be seen that the shear rate is dominating 
the fiber re-orientation compared to the stretch flow. Therefore, the shear- 
flow dominated flow in the commercial process simulation tool does not have 
the capability to predict the correct fiber orientations for the plug-flow dom- 
inated SMC flow. 
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Figure 4-6: Progression of the Young’s modulus over the plate length in and perpendicular 


to the flow direction (data from Mrkonjic et al. [Mrk15]) 


4.2. Compression molding simulation by using the 
CEL approach within Abaqus 


Due to the inaccuracies in the flow field and in the fiber orientation prediction 
of the commercial compression molding software, anew method needs to be 
developed to model the typical SMC plug flow. Therefore, the Coupled-Eu- 
lerian-Lagrangian (CEL) method within Abaqus is used, as it can describe a 
material flow as well as Fluid-Structure-Interactions (cf. section 2.5). The fol- 
lowing subsections introduce the modelling approaches and the implementa- 
tion in the Abaqus subroutine environment. Furthermore, the main process 
setting effects and the interaction between the SMC material and the rein- 
forcements is discussed. 
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4.2.1. Development of a deformation-rate based material 
modelling for SMC 


To overcome the disadvantages ofthe commercial compression molding soft- 
ware, a separate modelling of the core region and the lubrication layer is cru- 
cial. The lubrication layer develops due to the temperature gradient between 
the hot surface and the cold SMC material and is containing no fibers. Due to 
the relatively low thermal conductivity, this lubrication layer is only a few mi- 
crometers thick, depending on the exact process conditions [Lee87], [Sho14] 
(see Figure 3-22). Due to the relatively low viscosity in this lubrication layer 
compared with the viscosity in the core region, a high velocity gradient occurs 
and therefore a high shear deformation. Due to this effect, the core region 
shows a plug flow behavior, since no velocity gradient over the thickness oc- 
curs. These two flow behaviors cannot be described with the same material 
model and require separate formulations. Consequently, the stress in an ele- 
ment that comprises both regions can be expressed as 


o= -pl + Lrric + Orneo > (44) 


where p represents the pressure due to compression, 2¢,;, the friction stress 
in the lubrication layer and Grheoo the rheological stress in the core region. 
The modelling of these stresses is described in the following Subsections. 


Core region 


As shown in Section 3.3 and 3.8, the SMC material is a compressible Non- 
Newtonian fluid. Therefore, the Trouton ratio (see eq. (41)) between the 
shear and extensional viscosity is not valid. Due to this fact, the shear Nsnear 
and extensional viscosity Neyt need to be calculated separately. For the ex- 
tensional viscosity, the same basic formulation is used as in the rheological 
characterization modelling, and previous experiments show a good correla- 
tion (cf. Section 3.4): 
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IDIN” 
MDE Nps (Bt), (45) 
where Nps represents the plain strain extensional viscosity (comp. Section 2.2 
and 3.4). This equation needs to be extended to consider the temperature 
dependencies. Le Core et al. [Le 02] have examined the temperature depend- 
ence for the SMC paste for the extensional viscosity and achieved a good 
agreement with an Arrhenius equation. Therefore, the extensional power law 
function is extended with an Arrhenius function, so it can be written as 


nexe(D,T) = nps (Zt) exp (-» (+- 2). (46) 


To 
In this equation, Ty and b are fitted parameters. 


The extensional viscosity as well as the shear viscosity show an increase in 
their amplitude while the stress rate is reduced. This leads to quite large vis- 
cosities for infinitesimal low extensional and shear rates, which could lead to 
unphysical deformations. Therefore, a maximal extensional and shear viscos- 
ity is defined by using a minimal extension rate and a minimal shear rate, re- 
spectively. This behavior is observed for many materials and therefore, it is 
considered for example in the Cross-WLF viscosity model. These two mini- 
mum rates are considered as material parameters in the material card. 


For the shear viscosity, the same model is used as in the commercial com- 
pression molding tool (see eq. (42)). In contrast to the commercial process 
simulation tool, this model is acceptable, since the shear flow and, therefore, 
the shear viscosity has a minor influence on the total stress compared with 
the extensional based stresses. 


Based on this equation, the rheological stress can be written as: 


Orneo (D, T) = Nelong (D, T) dev(D) + shear (a, T, ý) Y, (47) 


88 


4.2. Compression molding simulation by using the CEL approach within Abaqus 


were dev(D) represents the deviatoric components of the deformation rate 
tensor D. 


To solve this rheological stress, the deformation rate tensor D needs to be 
calculated. Therefore, the deformation gradient F and the stretch tensor U 
are used to do the decomposition of the deformation rate tensor D. To do so, 
both tensors are used to calculate the velocity gradient L, the deformation 
rate D and the asymmetric part of the decomposition (spin matrix) W in the 
element and the global coordinate system. The main decomposition of the 
deformation gradient tensor F is given Figure 4-7. The deformation gradient 
is defined by deriving the coordinates of the current configuration xi from the 
coordinates of the reference configuration xo. In the case of the CEL method, 
the reference configuration (see Figure 2-10) is the configuration at the be- 
ginning of the increment. Therefore, the deformation gradient is defined as 


_ ak) 
Fij = ax,’ (48) 
where 7; (È, t) is the current position of a material point and X; the position 
of the considered material point in the undeformed output or the reference 
configuration. Due to the CEL formulation with the mapping step, this is equal 
to the configuration at the beginning of the increment. 
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Undeformed configuration 


configuration 
ez, 
Va 
e 
Figure 4-7: Polar decomposition of deformation gradient tensor F 
Based on this, the velocity gradient L, the deformation rate in the second in- 
termediate configuration D and the spin matrix W can be calculated by using 


the intermediate stage of the strain rate tensor &, the elastic part of the body 
rotation matrix R® and the symmetric plastic parts of the velocity gradient 


b= 22 pra (49) 


The spin matrix W and the intermediate stage of the strain rate tensor & are 
the decomposition of the velocity gradient L, defined as L = €+ W, and can 
therefore be expressed as: 


w=: (l-E) =i e-m. 60) 
(urn. ED 
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By dividing the deformation gradient tensor F and its composition into an 
elastic and a plastic part, the following correlations can be used: 


Re! = Fe! 3 uel”, (52) 
—pl 

D = R è- Re”, (53) 
D= R° D”. Re". (54) 


Beside this incompressible formulation, a formulation based on the determi- 
nant of the deformation gradient tensor is used, which represents the volu- 
metric change of an element. For that purpose, the compressible velocity gra- 
dient can be expressed as 


L 


Liso = Jamie 


(55) 
Based on this equation, the compressible second intermediate configuration 
of the deformation tensor D are calculated in the same way. This formulation 
of the deformation forms the basis to predict the material stresses and the 
fiber orientation evolution. 


Lubrication layer 


Due to the high thickness ratio between the lubrication layer and the core 
region in the order of three magnitudes, the thin layer should not be modeled 
as a separate phase. This would lead to a discretization which current compu- 
tation cluster cannot handle. Therefore, the stresses from the lubrication 
layer Zfric are modeled as a friction term at the mold surface. For that pur- 
pose, the same friction model as for the rheological characterization devel- 
oped by Abrams [Abr03] (comp. Section 3.4) is used: 


ee (ae a (56) 


vo vll all’ 
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where A represents the hydrodynamic friction coefficient, Ù the velocity ten- 
sor, m the power-law coefficient for the friction stress and ñ the surface nor- 
mal of the mold surface, which is perpendicular to the velocity tensor Ù. 


Both material formulations are defined as scalar functions. Therefore, the 
stresses need to be calculated in every element in the local coordinate system 
defined by the eigenvectors of the deformation rate tensor D and the velocity 
field Č, respectively. 


Compression pressure 


The pressure due compression p is calculated based on the incompressible 
hyperelastic specific deformation energy wọ. The general stress definition 
based on the three-dimensional continuum consideration is defined as: 


rn (57) 


eee dE 


with the Green-Lagrangian strain tensor E. To ensure a hyperelastic incom- 
pressible material formulation, the specific deformation energy Wg is defined 
as: 


Wo = ~~ (det(F) =i. (58) 


By deriving the specific deformation energy wọ and using the definition of the 
bulk modulus K of 


-y2 p? 
K=-v2= p2, (59) 


the compression pressure can be expressed as by using the resulting stress o 


p = ZIn(det(F))b*! , (60) 
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where b is the Left-Cauchy-Green tensor, which is defined as 


b=F-F. (61) 


4.2.2. Implementation of the SMC-Material model 


The material models are implemented within a user material subroutine 
(VUMAT), since the usage of a formulation based on a user equation of state 
(VUEOS) and a user viscosity subroutine (VUVISCOSITY) using the Trouton ra- 
tio. As described before, this fixed relation between shear and extensional 
viscosity needs to be extended for the SMC material formulation and could 
therefore not be used in this case. Furthermore, not all necessary input pa- 
rameters for the fiber orientation models are available with a combination of 
VUEOS and VUVISCOSITY, while they are available within a VUMAT. 


To implement the superposition of the three stress terms (cf. eq. (44)) most 
efficiently, the lubrication layer stress term should be implemented as a mold 
surface friction term. Since none of the standard friction models implemented 
in Abaqus Explicit has the correct form, a user subroutine friction model 
(VUFRICTION) needs to be implemented. Due to the Lagrangian and Eulerian 
step within the CEL method and their implementation within Abaqus, this 
subroutine is not supported within CEL. The same applies for the used inter- 
action subroutine (VUINTERACTION). Therefore, all stress terms have to be 
calculated in the VUMAT subroutine. 


The pressure stress and the rheological stress need to be calculated for all 
filled Euler elements, while the friction stress term is only active, if an SMC 
filled Eulerian element is in contact with a Lagrangian mold surface element. 
Within the VUMAT subroutine, this cannot be tested directly with build-in 
functions or queries. Therefore, the interaction information needs to be 
passed from other subroutines. As the usage of the VUFRICTION is not possi- 
ble, no other subroutine can pass this information directly to the VUMAT. To 
overcome this problem, a series of scripts and subroutines are implemented, 
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which identify the Eulerian elements that are in interaction by geometrical 
tests. Therefore, a self-implemented plug-in within Abaqus CAE initializes the 
geometrical information. This information is stored in individual state varia- 
bles (SDV) for every Eulerian element. Since the state variables are trans- 
ported with the fluid, latter are stored as a field variable in a second step, as 
they are not transported. Due to the definition of the initialization step, it is 
not possible to define the parameters directly as field variables for the 
EC3D8R element type in Abaqus CEL. This transfer from the SDV to the FLD is 
done for every new considered element within the user subroutine to rede- 
fine field variables (VUSDFLD). As this leads to a parallel processing of data in 
different subroutines and writing into the same array in different threads, the 
user-defined external subroutine VUEXTERNAL is used to secure the thread 
safety. These non-transported field variables are then used to compare the 
initialized information with the total motion of the moving mold side in the 
VUSDFLD. Therefore, the total motion of the moving mold side is passed in by 
the user defined amplitude subroutine (VUAMP). This information is then 
passed back to the VUMAT to activate the friction stresses for these Eulerian 
elements. 


A schematic illustration of the order and interactions of scripts (dark blue), 
subroutines (mint-colored) and internal information (black) is given in Figure 
4-8. Detailed information about single scripts and subroutines are given be- 
low. 
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Common model setup for CEL in Abaqus nitialization Plug-In 
7 Plug-In to define 


initial fiber 
Plug-In to initialize the geometrical BC orientations and 
and writing them into fields default SDVs 


BuIss39014-91d 


SDV init 


VUAMP 
Mold closing profile 


VUMAT 

Calculation of stresses 

* Core material 

* Lubrication layer 
Calculation of fiber 
orientation 
Calculation of stresses for 
reinforcements (optional) 


Speed vs. force control 


VUSDFLD 
Testing BCs 
Interaction of subroutines 


EL 


Const. geom. information 


VEXTERNALDB 
Thread safety 


Figure 4-8: Schematic illustration of the order and interaction of scripts (dark blue), sub- 
routines (green) and internal information (black) of the SMC process simulation 


Initialization Plug-In 


The lubrication layers occur on the moving and on the fixed mold surface. 
Therefore, different approaches are necessary to detect the Eulerian ele- 
ments, which are in interaction with the Lagrangian mold elements. Since the 
Eulerian elements have no relative velocity to the fixed mold surface, always 
the same elements are in interaction. Therefore, the flag as well as the inter- 
action area between the elements and the normal vector of the fixed mold 
surface are stored as state variables by the plug-in. For the interaction with 
the moving mold surface, the distance between the element and the moving 
mold surface in closing direction is necessary. When the nearest mold surface 
element in closing direction is detected, the minimum and the maximum dis- 
tance of this Eulerian element are stored in the state variables as well as the 
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normal vector of the corresponding element of the moving mold surface. To 
calculate the interaction area between the mold surface and the Eulerian el- 
ement in the VUMAT, the coordinates of the eight nodes of the Eulerian ele- 
ment as well as one reference point of the nearest mold surface are stored as 
state variables. To hand over the necessary input to the plug-in, a drop down 
menu is integrated in the user interface to choose the part instances for Euler 
mesh, moving mold and fixed mold (see Figure 4-9). In the next step, the clos- 
ing direction with the amplitude information is inserted for the detection of 
the moving mold element. 


To prevent a difference between the predefined closing direction and the one 
inserted in the user interface, the user-defined information is also used to 
create the boundary conditions. Since the amplitude of the moving mold 
boundary condition is controlled by a used amplitude subroutine (VUAMP), 
the initial amplitude needs to be defined. Therefore, a subroutine specific 
property value is used. Since this specific property cannot be assigned via 
Abaqus CAE, the keyword editor is used within this initialization plug-in. 


x 


This Plugin initalize the SMC VUMAT 


Specify closing direction: 


Choose Euler Part: | Modell: Euler hd 


Choose Moving Mold Part: Modell: mold 


Choose fixed Mold Part: Modell: mold-fixed 


write SDV and FLD Cancel 


Figure 4-9: Initialization Plug-in interface to initialize the geometrical information for the 
lubrication layer detection 


The second part of the plug-in can be used to initialize the fiber orientation of 
predefined element sets, which represent the initial charge. This can be done 
by assigning the fiber orientation to element sets or to each individual ele- 
ment by using the state variables. Since the initialization of the geometric in- 
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formation defines the state variables of each element, the initial fiber orien- 
tation distribution also needs to be defined element-wise. Since this defini- 
tion is not available within Abaqus CAE and in order to simplify the function 
for the user, the plug-in offers the definition of fiber orientation distribution 
for predefined element sets (see Figure 4-10). By using the keyword editor, 
the initial fiber orientation distribution is then written into the model. 


x 


Choose Model: ‘Model-Slice-D01_Rough_solid fl Define the components of the fiber orientation tensor> 


This Plugin initalize the SMC VUMAT J 


write FODs Cancel : |0 


| + x 
| Insert element set for the current FOD definition 


| elem_set: | ElemSet-1 


Figure 4-10: GUI to initialize the fiber orientation tensor for predefined element sets of the 
initial charge 


VUAMP - Closing speed control of the moving mold 


The user amplitude subroutine VUAMP is used to determine the closing speed 
of the current increment. The amplitude of the closing speed in the first in- 
crement is achieved from the Initialization Plug-in, by the property definition. 
For the further increments, the amplitude is calculated depending on the cho- 
sen closing mode. Since the closing direction is defined by the unity vector, 
the closing speed amplitude is independent of the closing direction. 


Besides the update of the amplitude, the VUAMP calculates the total displace- 
ment of the moving mold in closing direction from the beginning of the step. 
This information is needed in the VUMAT to identify the Eulerian elements, 
which are in interaction with the moving mold surface. This could not be done 
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directly within the user amplitude subroutine, as it is not accessing the single 
elements. Therefore, the current total movement is passed to the VUSDFLD 
to be processed thereby using a global array. Since the VUAMP has no direct 
output to state or field variables, the global array is the only possibility to 
transfer data. In this case, a global array is an array that can be freely assigned 
and called by all subroutines. 


Within the VUAMP, an updated amplitude for a constant deformation rate or 
a constant closing speed can be calculated. Up to now, no switch from speed 
to force control is implemented, but all interfaces for such a control are pre- 
defined within the subroutine. 


VUSDFLD - User defined state variables and field subroutine 


This subroutine has the two targets: to identify the Eulerian elements which 
are in interaction with the mold surface and to transfer the geometrical infor- 
mation from the state variables to the field variables. This is necessary to pre- 
vent the usage of averaged geometrical information, which could occur due 
to the advective transport and therefore mixture of state variables infor- 
mation. 


To identify the interaction elements, the current total displacement is 
achieved from the VUAMP and is compared with the two limits ofthe minimal 
and the maximal distance to the nearest mold surface point atthe initial state. 
By using a flag, this interaction information is passed to the VUMAT by the 
state variables. Since the VUSDFLD is called before the VUMAT, the transport 
of the state variables is not relevant. To prevent a simultaneous multiple us- 
age of the global array, the thread safety needs to be secured. Therefore, a 
mutex initialized within the user defined external program subroutine (VEX- 
TERNALDB) is used. This mutex blocks the global array for other processes, 
while one process is accessing and is writing input to it. 
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The geometrical information for the interaction needs to be stored as ele- 
ment-based information. Therefore, the state variables are used, as they are 
the only available element-based arrays in the initial step. The field variables 
can only be initialized at the nodes, which would lead to a mixture of infor- 
mation on the element level. Due to the effect that the state variables are 
transported by the advective terms, a mixture would occur leading to wrong 
geometrical information. Within the VUSDFLD, it is possible to define field 
variables at the integration points, which are the element centers for the Eu- 
lerian elements (EC3D8R or EC3D8RT). Therefore, all geometrical information 
is transferred from the state to the field variables during the first increment. 
Apart from that, since the CEL method is not considering every Eulerian ele- 
ment in every increment, it is checked, whether a new element is processed 
and, if necessary, the information is transferred. 


VUMAT - User material subroutine 


Within this VUMAT, the SMC material stresses and the fiber orientation dis- 
tributions are calculated. For the stress calculation, the superposition of the 
compressibility pressure p, the rheological stress O,ye. and the hydrody- 
namic friction stress &,,.;. is used (cf. eq. (44)). To calculate these stresses, the 
developed models from Section 4.2.1 are used. 


The compressibility pressure is based on the incompressible hyperelastic ma- 
terial formulation given in eq. (60). Therefore, the bulk modulus is imple- 
mented as a material property, which is stored in the material card. 


The rheological stress is calculated using an anisotropic viscosity tensor n de- 
pendent on the deformation gradient D, which is composed of the exten- 
sional and shear viscosity. Therefore, the extensional and the shear viscosity 
are using different models. For the shear viscosity, the reactive-viscosity 
model (see eq. 42) is used. To calculate the extensional viscosity, the temper- 
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ature dependent power-law function from eq. 46 is implemented. The coeffi- 
cients of the shear and extensional viscosity models are implemented as ma- 
terials properties, which can be easily modified within the material card in 
Abaqus CAE. 


The third stress component is the friction stress of the lubrication layer. As 
already described in the plug-in and VUSDFLD Section, this term is activated 
by a flag, which indicates, that the current Euler element is in interaction with 
the mold surface. In this case, the stress is calculated by using eq. 56. If the 
Euler element is in interaction with the fixed mold surface, the interaction 
area is precalculated and stored in the geometrical information in the field 
variables, before the simulation is started by the plug-in script. For the inter- 
action with the moving mold, the interaction has to be calculated during every 
increment. Therefore, the coordinates of the corresponding element nodes 
as well as a reference point of the surface and its surface normal are used. 
This information is also part of the element specific field variables. The two 
parameter of this model are also part of the material card within Abaqus CAE. 


Beside the stress update, the fiber orientation changes due to the flow kinet- 
ics are calculated. Since the fiber orientation is transported with the material, 
the advective transport properties of the state variables are used to store the 
components of the fiber orientation tensor and to ensure the correct 
transport. Due to this transport, the stochastically based constraint of 
tr(A) = 1 needs to be checked and restored at the beginning of every incre- 
ment, if necessary. 


After transferring the fiber orientation tensor into its eigenvector coordinate 
system, the fiber orientation change is calculated by using the Folgar-Tucker 
ARD-RSC model (cf. Section 2.4). To approximate the fourth order fiber orien- 
tation tensor A from the second order tensor A, the Invariant Based Optimal 
Fitting (IBOF) closure approximation is used (cf. Section 2.4). Therefore, the 
new fiber orientation tensor can be calculated by 


Anew = Aoa + Adt . (62) 
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After transferring this new fiber orientation into the global coordinate sys- 
tem, the tensor is stored in the state variables, which are passed to the next 
increment. 


4.2.3. Effects of simulation settings on numerical stability 
and accuracy 


Two parameters have a major influence on the simulation quality: the bulk 
viscosity and the time increment range. 


The bulk viscosity introduces an extra damping to the simulation, associated 
with the volumetric strain. Although it is intended for high-speed events, it 
also has an influence on the Fluid-Structure-Interaction in the SMC simula- 
tion. Especially in the contact regions between the tool surface and the SMC 
material, a high influence of the bulk viscosity on the penetration can be ob- 
served. Due to this penetration, a parameter study based on the default val- 
ues is made. By increasing both, the linear and the quadratic bulk viscosity 
parameter, a better interaction surface with a lower penetration can be 
achieved. Three results with the default, an increased and a decreased linear 
and quadratic bulk viscosity parameter set by one decade is plotted in Figure 
4-11. Furthermore, the resulting total pressure is given in this plot. For the 
reduced and the original linear and quadratic bulk viscosity parameter set, 
the pressure is nearly the same. For the increased parameters, the pressure 
is significantly higher and the simulation is terminated due to numerical prob- 
lems. Due to these effects, an increase above the factor of 10 of the linear 
and the quadratic bulk viscosity parameter reduces the penetration, but leads 
to numerical instabilities in the simulation and should therefore be avoided. 
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S_ASSEMBLY_EULER-1_SMC_VUMAT-1, Pressure 


S_ASSEMBLY_EULER-1_SMC_VUMAT-1, Pressure 


JBLY_EULER-1_SMC_VUMAT-1, Prexsute 


Figure 4-11: Evaluation of the penetration of the SMC through the mold surface depending 
on the bulk viscosity parameter: reduced (top), original (middle) and increased 
parameter set (bottom) 


The second important setting for a Coupled-Eulerian-Lagrangian simulation is 
the range of the valid time increment range. Typically, this increment value is 
calculated automatically by using the element maximum eigenvalue and the 
speed of sound in the SMC material. Depending on the complexity of the cav- 
ity and the mesh quality, this increment can be too large. Therefore, a maxi- 
mum increment is recommended in the range of 1e-4 to 1e-5 s to ensure the 
correct considering of local deformation effects. 
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4.2.4. Consideration of local reinforcements 


The possibility to consider local reinforcements is one of the advantages of 
the CEL compression molding approach. Since the CEL formulation divides the 
increment into a Lagrangian and an Eulerian step (cf. Section 2.5), the inter- 
action of the SMC and a local reinforcement is considered in the Lagrangian 
step. In this step, the Fluid-Structure-Interaction is simplified to a Structure- 
Structure-Interaction using an extrapolated surface for the SMC material (cf. 
Section 2.5). Therefore, the tangential and the normal contact behavior need 
to be defined. Furthermore, the thermal conductivity could be defined, ifthe 
local reinforcement and the SMC are placed into the mold with different tem- 
peratures. Since the basic contact model used for the interaction between the 
SMC material and the tools is considering the interaction of every section el- 
ement with every other section, the reinforcement-SMC and reinforcement- 
tool-interaction need to be defined by assigning contact surface pairs. 


Since the material is flowing around the reinforcements and due to the pen- 
etration method for the interaction, it is recommended to use solid elements 
for the reinforcements. This ensures the distinctive determination of the sur- 
faces and their normal, as well as the flow directions at which they are lo- 
cated. Otherwise, the direction of the restoring force is not necessarily cor- 
rectly defined (see Figure 2-13). 


In most cases, the local reinforcement reduces the gap height significant. 
Therefore, a special focus on the mesh near the reinforcement is necessary. 
To predict a correct flow, at least four elements over the thickness are re- 
quired. This leads to an essential mesh refinement around the reinforcement 
and its probable moving direction to ensure the minimum element numbers 
during the whole molding process. 


For the material properties of the reinforcements, every available formula- 
tion defined for solid elements can be used. Furthermore, self-written user 
material subroutines (VUMAT) can be used. If this is done, a survey of the 
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material sections within both VUMATs is required to clearly identify the rela- 
tion between the element and the corresponding material. Such a VUMAT 
can be used to consider the material behavior of unidirectional tape rein- 
forcement and their fiber reorientation and fiber volume change. 


After the individual property assignment, the interaction properties need to 
be defined. Compared to the standard SMC compression molding simulation, 
interactions between the SMC and the local reinforcements and if necessary 
between the reinforcement and the molding tool have to be modelled. Exam- 
ples of such a consideration of isotropic local reinforcements by using a stand- 
ard elastic material card are presented and discussed in Section 4.3.2. 


4.3. Validation of hybrid SMC process simulation 


To validate the SMC model developed and implemented in Sections 4.2.1 and 
4.2.2, two different parts are used. On the one hand, the pressure distribution 
is compared with the results from the characterization experiments from Sec- 
tion 3.2. Beside the pressure, the fiber orientation distribution of the CT-scans 
(Section 4.1) is used as a validation criterion. On the other hand, the IRTG 
reference structure with local unidirectional reinforcements is simulated to 
evaluate the flow behavior of the SMC. For this part, the fiber orientation and 
the reinforcement displacement are compared with the experiments. 


4.3.1. 2D flow within the rheological tool 


4.3.1.1. Model configuration 


In this first validation step, the rheological tool presented in Section 3.1 and 
the pressure results shown in 3.2 are used. Therefore, the local effects and 
pressure changes should be resolved within the CEL model of the tool. In or- 
der to stay with a reasonable number of elements, only a section of the rhe- 
ological tool is simulated. To ensure the same flow as in the real process, sym- 
metrical boundary conditions are assumed. This is permissible since the 


104 


4.3. Validation of hybrid SMC process simulation 


pressure sensors are aligned in the middle of the tool. The model has a width 
of 10 mm, which is resolved with three elements. Three elements are used as 
a compromise to avoid fully constraint elements, if one element is used, on 
the one side, and to keep the number of elements low on the other side. In 
flow direction, the element edge length is set to 0.5 mm to precisely describe 
the flow front. This leads to a total number of 240 000 Eulerian elements. 


Figure 4-12: Boundary conditions of the Lagrangian Elements for the rheological tool CEL 
model with the initial charge (beige) configuration of the semi-structural SMCs 


For the validation, the Class-A SMC, the Vinylester SMC and the UPPH carbon 
fiber SMC are simulated. The Low-Density Class-A SMC is omitted, since it has 
nearly the same material properties as the standard Class-A SMC (see Table 
3-6). 


For the process settings, the same values as for the characterization methods 
are used. To reduce the amount of simulations, only one closing speed of 
h = 1 mm s” is considered. To ensure the correct viscosity of the SMC mate- 
rial, the initial temperature is set as 303.15 K (30 °C), which is also the refer- 
ence temperature in the Arrhenius equation (see eg. (46)). 


4.3.1.2. Simulation results 


To validate the simulation models with the experiments, various simulations 
were carried out. Therefore, the different initial charge configurations for the 
different material formulations were considered as well as different mold 
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closing configurations with constant closing speeds and constant deformation 
rates. 


During the simulation runs, an overall problem with the rheological tool 
model has emerged. Due to the boundary conditions with the moving tool 
side and the fixed tool side, only two degrees of freedom are left. These are 
the flow in x-direction and the compressibility of the SMC. Depending on the 
stress magnitudes of the rheological stress due to the deformation rate (see 
eq. 47) and the compression pressure (see eq. 60) this leads to a high com- 
pression of the material. Based on this compression and the resulting stored 
energy, little numerical variation within the model leads to high local stresses 
when the volume is expanding again, and therefore, to high local defor- 
mations. In some cases, this leads to very small time increments which are 
too small to simulate the whole compression process within an acceptable 
computational time. In other cases, these deformations lead to too high de- 
formations and an error in the simulation. To identify the reason for these 
issues, a mesh refinement was done. This study showed no improvements 
regarding these local effects, so that mesh effects could be ruled out. In a next 
step, the rheological tool model was tested by using the same models as for 
the Moldflow simulations as well as standard implemented models from 
Abaqus. To utilize the models used by Moldflow, the Castro-Macosko model 
was implemented in a viscosity user subroutine (VUVISCOSITY) and the 2-Do- 
main Tait Model in an user defined equation of state subroutine (VUEOS). 
Within these subroutines, the same material parameters as in the Moldflow 
simulations were used. This simulation also shows very small time increments 
in the same range as the SMC specific subroutines but slightly smaller. Fur- 
thermore, also this simulation gives an error after 5e-5 s due to local density 
changes (see Figure 4-13 and the corresponding high stresses, which lead to 
too high element deformations. 
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Figure 4-13: Local density change within the Castro-Macosko model with the 2-Domain-Tait 
model 


Additionally the standard Abaqus Mie-Grüneisen [Bae] equations of state and 
a constant viscosity in the same range as before were used. This simulation 
has a six decades smaller time increment, which means, that no results could 
be generated. The main reason for the very small time increments can be 
traced back to the viscosity value. By reducing the viscosity by one decade, 
the time increment is also reduced by one decade. This becomes obvious 
when looking at eq. 47, were the viscosity has a linear correlation with the 
stress and therefore with the deformation. 


By comparing the results of the SMC specific material model, the Castro-Ma- 
cosko model combined the 2-Domain-Tait model and the standard Mie-Grü- 
neisen model with a constant viscosity, the advantages of the new developed 
SMC specific material model becomes clear. Due to the high viscosities of the 
SMC during the compression molding process small time increments cannot 
be avoided. Due to the separation of the extensional and shear viscosity, the 
SMC specific material model can handle them better than the standard 
models. 
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Since previous tests showed the capability ofthe presented SMC compression 
molding approach by using the CEL method for 3D flows, this is only a problem 
of this restricted model. One of these test results is given in Figure 4-14. 


Figure 4-14: Results of a test mold filling simulation during the development process 


To find the reason for the issues, the SMC specific material stresses (rheolog- 
ical and friction stress) and the compression pressure of a simulation run with 
a VE-SMC material at the sensor position 1 and 2 are shown in Figure 4-15. In 
this Figure, it can be recognized that the first typical stresses emerge after 
0.075 s step time. This is due to the initial gap between the SMC and moving 
mold side of 0.05 mm to avoid initial problems. Therefore, the moving mold 
is in contact with the material after 0.05 s and the first material flow occurs 
after this. After the flow begins, the compression pressure reaches a 20 time 
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higher level than the rheological stresses. After approx. 0.7 s step time the 
first inaccuracies in the simulation occurs and the stress values at Sensor 1 
are dropping. These phenomena are occurring all along the model and cannot 
be traced back to geometrical or numerical reasons. 
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Figure 4-15: Rheological, friction and compression pressure stress for a VE-SMC material 


formulation at sensor position 1and 2 


Having a detailed look at the rheological and friction stress a slide oscillation 
of the rheological stresses can be observed. This is due to the change of the 
deformation rate tensor D due to the compression changes. It can be also 
observed, that the friction stresses from the simulation is different from the 
analytical model. The friction stresses at Sensor 1 is zero over the given step 
time. This difference occurs, since the analytical model uses an integral meas- 
ure over the complete flow length, while within the simulation discrete 
stresses are calculated. Within the simulation model, the total pressure is the 
integral equivalent to the superposition of the rheological and the friction 
stress in the analytical model. Here, the discrete friction stress at Sensor 1 is 
lower than the stress at Sensor 2, since no axial flow occurs at the end of the 
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rheological tool, the friction stress needs to be zero (see eq. 56). For the fric- 
tion stress at Sensor 2 a slide increase can be observed, since the velocity of 
the material is increasing at a fixed point during the compression molding 
process. 


Nevertheless, the simulation results show the typical plug-flow. Three short- 
shoots are shown in Figure 4-16. Within this Figure, a slide ahead running of 
the SMC material near the mold surface can be observed. This is a typical SMC 
specific material behavior, which was already observed by Barone and Caulk 
[Bar85] (see Figure 2-4). 
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Figure 4-16: Evaluation of the flow front with the compression stresses at t=0 s (top), 
t = 0.375 s (middle) and t = 0.75 s (bottom) 


Since the simulations are not finishing, no validation with the experiments can 
be performed. Nevertheless, to show the capability of the new developed 
material model, the different material formulations are compared with each 
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other. Therefore, a quantitative evaluation can be performed, since flow ef- 
fects due to the different initial conditions are excluded. To do so, the 
Class-A configuration with an initial charge length of lo = 300 mm is used. 


These simulations have the issues with the compression behavior like the pre- 
vious presented VE-SMC simulation results. Nevertheless, a sensitivity to the 
material parameter can be observed, when comparing the stable step times. 
The Class-A SMC has the shortest stable step time, while the VE-SMC runs 
twice as long and the UPPH-SMC thrice the time. This result, together with 
the pre-studies, shows a special sensitivity to the rheological stress power- 
law coefficient n. Since this is higher for the Class-A SMC compared to the 
semi-structural SMCs (see Table 3-6), a higher sensitivity to changes of the 
deformation rate tensor D occurs. This can also be seen in the lower stable 
time increment during the simulation. Since the compression leads to such 
local changes, it is to be expected that within the Class-A SMC simulation high 
element distortions occur earlier. 


To compare the simulation results, the stress distributions for the Class-A 
SMC is given in Figure 4-17, for the VE SMC in Figure 4-18 and for the carbon 
fiber UPPH SMC in Figure 4-19. All these three pressure curves show the SMC 
typical pattern with the high increase at the beginning due to the rheological 
stress and the threshold after this first increase. The semi-structural SMC for- 
mulations show some overshoots at the beginning, but converge to the 
threshold after a short time. This is once again an effect of the compression 
stress and the resulting local deformation rate change. 
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Figure 4-17: Simulated stress distribution for Sensor 1 to 3 for the Class-A-SMC material 
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Figure 4-18: Simulated stress distribution for Sensor 1 to 3 for the VE-SMC material 
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Figure 4-19: Simulated stress distribution for Sensor 1 to 3 for the carbon UPPH-SMC 
material 


Comparing these three figures it can be observed, that the correct tendencies 
of the pressure peaks compared to the material parameters (see Table 3-6) 
are obtained. The Class-A SMC has a higher extensional viscosity and a higher 
rheological power-law coefficient than the UPPH carbon fiber SMC. Conse- 
quently, the predicted stress is also higher. The VE-SMC has the same material 
parameter as the UPPH carbon fiber SMC, but a higher extensional viscosity. 
Therefore, the predicted stress is higher. 


Nevertheless, comparing the VE-SMC and the carbon UPPH SMC, gives the 
reasons for the difference to the experiments. Comparing the experimental 
pressure curves for the VE-SMC (see Figure 3-5) and the UPPH SMC (see Fig- 
ure 3-7) and on the extensional viscosities given in Table 3-6, the VE-SMC 
should have a three times higher stress threshold than the carbon UPPH SMC 
while having the same deformation rates. Nevertheless, the simulation shows 
has a 30 % higher pressure for the VE-SMC formulation. Having a detailed look 
at the simulation results, Figure 4-20 shows a comparison of the compression 
based on the determinant of the deformation gradient (det(F)) and Figure 
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4-21 a comparison of the displacement in flow direction for the two semi- 
structural SMC materials. These two figures show the effect of the compres- 
sion on the displacements and therefore also on the deformation rate (see 
eq. 49 to 54). Already after 0.025 s simulation time, the VE SMC process sim- 
ulation shows a higher compression (lower determinant of the deformation 
gradient) compared to the UPPH SMC process simulation. This can be at- 
tributed to the higher rheological stresses due to the higher reference plain 
strain extensional viscosity of the VE SMC. Since the compression leads to a 
reduction of the volume, some material is flowing opposed to the flow direc- 
tion near the end of the cavity (negative material displacement at the left end 
in Figure 4-20 and Figure 4-21). This change of the idealized constant defor- 
mation along the material leads to a reduced deformation rate tensor D and 
therefore to a reduction of the local rheological stresses. Furthermore, an os- 
cillating displacement near the flow front can be observed. In this area, the 
displacement is changing between zero and the maximum value. Since these 
high gradients are not occurring once, but over the half material, this is once 
again an effect of minimal numerical changes in the compression at an incre- 
ment. 
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Figure 4-20: Comparison of the compressibility (det(F), see eq. 60) of the carbon UPPH SMC 
(top) and VE SMC (bottom) after 0.025 s simulation time 
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Figure 4-21: Comparison of the displacement in flow direction (U1) of the carbon UPPH SMC 
(top) and VE SMC (bottom) after 0.025 s simulation time 


Since the presented simulation results are made with the Class-A SMC initial 
configuration, a comparison of the VE-SMC simulation with the experiments 
can be done roughly. Therefore, both pressure curves are combined in Figure 
4-22. The experimental pressure curves have the compaction phase at the 
beginning, were the entrapped air is released and inhomogeneities in the 
thickness are smoothed. Since the simulation starts with ideal conditions, this 
is not modeled and therefore different times are needed until the rheological 
stress threshold is reached. If this is considered within the definition of the 
Initial Charge thickness, this has no effect for the further simulation. Further- 
more, the thickness inhomogeneities are the reason why the experimental 
curves are not on top of each other in this phase. Nevertheless, this Figure 
shows the correct stress prediction of the rheological stress during this phase. 
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Figure 4-22: Comparison of the simulation results with the experiments of the Class-A ma- 


terial formulation 


These process simulation results show the possibilities of the new developed 
SMC process simulation based on the CEL approach. Comparing the simula- 
tion results from the commercial compression molding simulation with the 
SMC-specific material model results within CELthe advantages of the new ap- 
proach becomes obvious. Comparing the flow fronts of the commercial sim- 
ulation tool (see Figure 4-5) with the SMC-specific material model (see Figure 
4-21) andthe pressure curves (see Figure 4-2 and Figure 4-15) the correct flow 
kinetic with the SMC specific plug-flow can be predicted with the new ap- 
proach. Due to this flow kinetics with the elongation based core region and 
the resin-rich lubrication layer, the pressure distribution can be predicted bet- 
ter due to the superposition o the rheological and friction stresses. Further- 
more, the reliable flow prediction leads to a better prediction of the Fiber 
Orientation Distribution, since they only depends on the flow fields. 


Nevertheless, due to the issues with the compression behavior during the 
compression molding simulation, these simulations were not able to finish 
without an error. Therefore, better suitable compression pressure models 
needs to be developed, which can fulfill the criteria for the Coupled-Eulerian- 
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Lagrangian approach and the SMC specific compression behavior. Therefore, 
some extensions, which might replace the scalar bulk modulus by a tensor, 
could help to overcome the current issues. Such a tensor should not lead to 
an anisotropic bulk modulus but should describe the bulk modulus as a func- 
tion of stresses, compression state and other state variables. 


4.3.2. IRTG reference structure without ribs 


4.3.2.1. Model configuration 


Within the International Research Training Group (IRTG) GRK2078, a refer- 
ence structure was developed to show the different effects of the local unidi- 
rectional reinforcements in combination with SMC. This reference structure 
applies the rheological tool from Section 3.1, which has a replaceable insert 
in the middle section. By using this insert, a beading is introduced into the flat 
plate (see Figure 4-23). The beading angles are chosen differently in order to 
achieve asymmetric deformation of the molded part in the basic configura- 
tion in a four-point bending test. The beading angles towards the long plate 
ends are defined as 45°, while the other angles are set to 60° (positive y-di- 
rection) and 40° (negative y-direction). To create more homogeneous ten- 
sions, two unidirectional reinforcements with slightly different widths were 
inserted. The shape and the position of these two unidirectional reinforce- 
ments are also given in Figure 4-23 (green part). 
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Figure 4-23: IRTG reference structure with the position of the local reinforcements (green) 


The preformed initial charge was placed over the beadings. The top layer of 
the initial charge had a cutout for the unidirectional reinforcements (see Fig- 
ure 4-24), so they are fixated at the beginning of the molding process (see 
Figure 4-24). The mold closing speed was set to 1 mm s“ and the maximum 
press pressure to 2500 KN. Due to the inserts, the pressure measuring system 
used for the rheological characterization could not be used during the part 
manufacturing. 


In the real part manufacturing, the semi-structural UPPH carbon and glass fi- 
ber SMC (cf. Section 3.2.4 and 3.2.5) were used as SMC materials. Due to the 
different and uncommon flow behavior of the UPPH glass fiber SMC (cf. Sec- 
tion 3.7), only the carbon fiber SMC is simulated. For the process settings, a 
closing speed of h = 1 mm s* were used. The initial charge had a temperature 
of 303.15 K (30 °C), while the tool was tempered to 418.15 K (145 °C). 
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Figure 4-24: Initial configuration of the process simulation of the IRTG reference structure 
with the local unidirectional reinforcements (green) and the initial charge 
(beige) 


To achieve a mesh with Eulerian elements, which are mainly aligned in flow 
direction and which have a high mesh quality, a block-structured meshing 
method is used. Therefore, the Eulerian domain is predefined within a CAD 
model by defining the maximum mold opening distance and the necessary 
element overlap over the Lagrangian mold surfaces (cf. Section 2.5). To en- 
sure a high mesh resolution in the area of the beading, while having a reduced 
resolution in the flat area, an O-Grid is inserted (see Figure 4-25) within ICEM 
CFD. Furthermore, a mesh refinement is introduced in the thickness direction 
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to ensure the minimum element number at the end of the molding simula- 
tion. Therefore, different layers with different element heights and asmooth- 
ing via biasing is introduced. After exporting and convertingthe mesh, it could 
be used within Abaqus. 


Layered 
Thickness 


Figure 4-25: Blocking of the Eulerian mesh of the IRTG reference structure 


To model the precured local unidirectional reinforcements, homogenized ma- 
terial properties based on a Cantilever test (DIN 53362) provided by the Insti- 
tute of Production Science (wbk) were used. The interaction between the 
mold surface and the reinforcement is assumed to be rough, so that no slip 
will occur, if both are in contact. 


The initial charge of the SMC is formed to the 3D geometry of the beading. 
Therefore, the fiber orientation is not homogeneously quasi transversal iso- 
tropic in the global x-y plane. To predefine the correct fiber orientations in 
the flanks of the beading, a script is written to subdivide the element set pro- 
vided by the volume fraction tool, which defines the position of the initial 
charge within the Eulerian elements. After getting the sets for the initial 
charge elements in the different flank regions, the correct rotated quasi trans- 
versal isotropic fiber orientation tensor could be predefined in the model. 


Due to the high number of elements and the relatively long flow length, the 
calculation with the same material parameters and process condition as in 
the experiments would take more than one month. Therefore, the lubrication 
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layer is not considered in this model, as the effect for such a one-dimensional 
flow can be neglected, since no pressure controlled closing speed is applied 
and no pressure validation is done. Furthermore, a mass scaling is done, 
which decreases the computational time dramatically, but slightly increases 
the compressibility of the material at the beginning of the compression mold- 
ing process. According to numerical pre-studies, these simplifications effect 
the result in an acceptable range compared to the amount of saved compu- 
tational time. 


4.3.2.2. Simulation results 


To validate the simulation results quantitatively, the patch displacement is 
the most important criteria. Although no UD-patch material model was used, 
but an isotropic one, the propensity to patch displacement or deformation 
can be evaluated. A representative picture of the final plate with the patches 
is given in Figure 4-26. This picture shows a configuration with the UPPH glass 
fiber SMC with continuous carbon fiber patches, also based on the UPPH 
resin. The patches are made out of two layers with a [0, 90] configuration. By 
evaluating the displacements of the local patches of the six manufactured 
samples, almost no displacement was observed in the main area (see Figure 
4-26). Only at the end of the tapes, were the thin arms are, some displace- 
ment could be observed (see Figure 4-26 bottom). In some cases, a small ro- 
tation of the whole patch could be observed in the final part, which was 
traced back to little misalignments during the stacking of the tapes and the 
initial charge due to the manual process. Since these rotations are initially 
introduced, also no displacement occurs during the process. 


121 


4. Process simulation of hybrid SMC 


Figure 4-26: Photos of the molded reference structure with local reinforcements with no 
deformation of the tapes (top) and slide deformation at the end of the tapes 
(bottom) [But17] 
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To compare the simulation results with the experiments, the patch displace- 
ments divided into the three coordinate directions are given in Figure 4-27, 
Figure 4-28 and Figure 4-29. The in-plane directions (x and y), plotted in Figure 
4-27 and Figure 4-28, show nearly no displacements of the patches. Since the 
mold-opening gap in the flanks is changing differently than the planar areas, 
a forming of the tapes into the final configuration can be observed. At the end 
of the tapes, were they are very small, further displacements can be ob- 
served. For the lower tape, this end region is elongated in flow direction up 
to 1.35 mm and shows a small deflection of 0.5 mm. For these in-plane dis- 
placements it can be observed, that both tapes are behaving differently. This 
is due to the different shapes of the tapes as well as the different angles of 
the flanks. This effect was also observed during the experiments, since the 
non-symmetric geometry of the reference structure leads to anon-symmetric 
flow. The prediction of these different deformation behaviors of both tapes 


shows the advantages of the SMC specific material model with the Fluid- 
Structure-Interaction. 


Figure 4-27: Displacement of the tapes during the compression molding process in flow di- 
rection 
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Figure 4-28: Displacement of the tapes during the compression molding process in-plane 
perpendicular to the flow direction 


In Figure 4-29 the displacement of the tapes in closing direction is given. Since 
the tapes are placed on top of the initial charge, they need to follow the upper 
mold movement. Therefore, they need to push back the SMC material under 
it without being compressed. The simulation results show that the planar re- 
gions of the tapes are following the displacement, while in the flank regions 
higher displacements can be observed. This is once again a forming process 
since the initial tape shape is not exactly the finale one. This forming in the 
flank area leads to a higher deformation of the lower patch in the main in- 
plane area (in Figure 4-29 at the very bottom) due to the relative stiff patch 
material. As a consequence, this displacements leads to a higher deformation 
of the SMC material and therefore to a higher pressure under this patch area 
(see Figure 4-30). The prediction of such an interaction is a clear advantage 
towards all available SMC process simulation tools. 
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Figure 4-29: Displacement ofthe tapes during the compression molding process in compres- 
sion direction (bottom) 
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Figure 4-30: Pressure of the SMC under the lower patch due to the patch deformation 


Comparing the experiments with the process simulation results in terms of 
the patch displacement, good agreements can be achieved, since the tape is 
following the tool movement in closing direction, but no displacement in the 
in-plane direction occurs. Beside the displacements it can be also observed, 
that no SMC material is flowing between the tool side and the patches. This 
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indicates that the interaction properties and simulation setting (see Section 
4.2.3) are also valid for this type of simulation with further Lagrangian parts. 


4.4. Summary process simulation 


In this chapter, the inaccuracies of commercial compression molding software 
for SMC were shown. To overcome these problems, a new simulation method 
based on the Coupled-Eulerian-Lagrangian approach is developed and pre- 
sented. This approach has the capability to predict the stresses of the SMC- 
specific core region and lubrication layer. Furthermore, the interaction with 
local reinforcements or inserts is possible. During the validation, these ap- 
proaches and their functionality were generally shown. By comparing the first 
threshold level of the rheological stresses at the beginning of the compression 
molding process for the Class-A SMC formulation, the correct pressure could 
be predicted in the simulation. Due to the compression behavior of the simu- 
lated material, the match for the semi-structural SMC formulations were not 
as good as the one for the Class-A SMC formulations. In the second step, the 
interaction of local reinforcements with the SMC was shown. By using the 
IRTG reference structure the correct displacement in flow direction was pre- 
dicted as well as the forming process in the closing direction. 


These two examples show the advantages of the new developed SMC specific 
material model. The model has the capability to predict the correct SMC spe- 
cific pressure distribution with its rheological and friction parts. Furthermore, 
the full interaction of the SMC with local reinforcements was shown, which is 
currently a unique method for Sheet Molding Compound compression 
molding. 
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Due to the fiber reorientation during the compression molding process, the 
consideration of the final fiber orientation distribution is crucial for the struc- 
tural simulation. Since this also affects the design, a holistic approach is nec- 
essary. Therefore, the Institute of Vehicle System Technology (FAST) — Light- 
weight Technology is focusing on the CAE-Chain, which combines different 
simulation tools and enables the transfer of information between them. 
Therefore, the single FE simulation tools, which are used during the steps De- 
sign, Optimization, Process Simulation, Mapping and Structural Simulation, 
are connected. In this CAE-chain, each step has criteria that must be fulfilled 
before results are passed to the next step. Otherwise, a process internal op- 
timization tries to fulfill the criteria by changing free parameters. If this is not 
possible or if the internal optimization fails, the previous step is done again 
with a varied parameter set. [Hoh16b] 


The structure of the CAE chain developed at KIT-FAST is presented in Figure 
5-1. 
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Figure 5-1: Steps of the CAE-Chain at KIT-FAST [Hoh16b] 


To transfer the information from the short or long fiber process simulation 
(injection or compression molding) to the structural simulation, a plug-in 
called Fibermap is developed within Abaqus, which can export and process 
the fiber orientation distributions. Furthermore, the mapping, clustering and 
a homogenization are part of the plug-in. All these steps of the Fibermap plug- 
in can be executed via a GUI within Abaqus (see Figure 5-2) or via single scripts 
by using the Abaqus command line. 
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Plug-ins Help %? 


Toolboxes _ | Ja | Qy Liste Cali g 
— rE — 
1 Export Abaqus Model 
InitSMC-VUMAT > 2 Convert Moldflow Results 
Tools Dean 3 Mapping 
About Plug-ins... 4 Import and Cluster Mapped Results 
5 Generate Materials 
6 Export Results 
7 Export CEL Results 
Figure 5-2: GUI selection tree of the Fibermap plug-in within Abaqus 


The processing of the mesh and result information is done using the neutral 
transfer format VTK. The VTK-file uses the ASCI encoding and consists of the 
mesh with the corresponding nodes and elements as well as the result infor- 
mation. This information is assigned to attributes. Therefore, various differ- 
ent types of results can be stored in a single VTK-file. In the following sections, 
the main steps of the CAE-chain provided by the Fibermap plugin are de- 
scribed. Therefore, the steps export of results, mapping and homogenization 
are discussed and the used methods are explained. 


5.1. Export of results 


The processing of the fiber orientation is realized for the simulation with the 
commercial process simulation tool Moldflow as well as for the internally ex- 
tended Abaqus CEL process simulation. Due to the different input file formats, 
two different scrips are necessary. The subscript “Convert Moldflow Results” 
uses the exported Patran mesh and the element based fiber orientation as 
input (see Figure 5-3). These two files have to be generated within the com- 
mercial process simulation tool before using the script. In some cases, not- 
filled elements can cause some issues for the post processing. Therefore, the 
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“Autoremove” function is inserted to give the user the option to remove 
these unfilled elements. As the outcome of the script, a VTK file including the 
fiber orientation distribution is generated. In this VTK file, the fiber orienta- 
tion distribution is defined in the tensor notation. Therefore, a 3x 3 result 
matrix is given for every element. 


This Plugin converts Moldflow results to VTK 


Geometry (*pat): 


Orientations (*ele): | 


Autoremove: True ir) 


Convert to VTK Cancel 


Figure 5-3: GUI of the VTK converter for the process simulation with Moldflow 


For the export of the fiber orientation from the process simulation with 
Abaqus CEL, the VTK file is generated from the result (odb) file, which needs 
to be opened in Abaqus CAE. Since the fiber orientation distribution is stored 
as a State variable (cf. Section 4.2.2), the Eulerian part instance and the Eu- 
lerian volume fraction (EVF) Section need to be chosen. These inputs enable 
an internal access to the corresponding state variables within Abaqus CAE. 
Due to the parallel computing, the element numbering provided by Abaqus is 
not in accordance with the VTK convention. Therefore, a renumbering of the 
elements is done to ensure a valid output. The generated VTK file is also using 
the tensor notation of the fiber orientation results. 


To enable an easy generation of the VTK file from the CEL process simulation 
results, the Export CEL Results GUI from the Fibermap Plug-in can be used. 
This GUI provides all available jobs, steps, instances and Eulerian Sections 
within separate dropdown menus to avoid wrong input to the export script 
(see Figure 5-4). 
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m 
= Export CEL Results to IRTG compliant VTK 8 


This Plugin exports an Abaqus Part with FOD results to IRTG compliant VTK. 


Clustering Result (*_labeled.csv): | C:\Temp\log_labeled.c| [Select] 


4 


Job and step: | E:/1513584254_Rheo_tool_arl_noFric_RMP_Mxe2/Rheo_tool_arl_noFric_RMP_Mxe2.odb+ Step-1 


Choose instance: | EULER-1 hd 


SDV: | ASSEMBLY_EULER-1_SMC_VUMAT-1 |x 


Figure 5-4: GUI of the VTK converter for the process simulation with Abaqus CEL 


Beside the processing of Abaqus CEL results, it is also possible to convert re- 
sults from Moldflow. Therefore, the patran file and the element based result 
information need to be exported from Moldflow. These two files can then be 
converted to a VTK-file within the Fibermap Plugin. 


5.2. Mapping of process simulation results 


In the most cases, the structural simulation uses a different mesh than the 
process simulation. Therefore, the information needs to be mapped from the 
one to the other mesh. This is done by mapping tools. Within the CAE chain, 
the mapping is done by using MapLib, which is part of the MpCCl Multi-Phys- 
ics Software Solutions provided by the Fraunhofer Institute for Algorithms 
and Scientific Computing (SCAI). This script-based tool needs two VTK files, 
the source mesh with the result information and the target mesh. While the 
VTK file with the result information is generated by the export results scripts, 
the target mesh VTK file needs to be generated separately. If the target mesh 
is available within Abaqus, the Fibermap Plug-in offers an exporting tool (see 
Figure 5-5). By choosing the desired part in an opened model, the VTK file of 
the target mesh is generated, with the correct element and node order con- 
vention. 
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x 


| This Plugin exports the geometry of an Abaqus Part. 


Choose Part: |Model-1: zugprobe e] 


Export to VTK Cancel 


Figure 5-5: GUI of the export script for the structural simulation mesh 


Beside the two VTK meshes, the MapLib mapper needs a configuration input. 
Within this configuration file, the mapping algorithms as well as the neighbor- 
hood search options are defined. These options have a significant effect on 
the mapping quality and need to be chosen depending on the structural mesh 
type (solid or shell), see below. All inputs can be passed to the mapping tool 
by the Mapping GUI (see Figure 5-6). As a result, a new VTK file based on the 
structural target mesh input with the mapped information from the source 
mesh is generated. 


This Plugin maps Moldflow results to an Abaqus Mesh 


mr m 
Target | [Select vrK 


Configuration: c:\Users\uadsy\Desktoy Select config.rc 


Figure 5-6: GUI for the input of the MpCCI MapLib mapping tool 


To identify the best mapping algorithm, a generic study with finer and 
rougher target meshes is performed. An exemplary mapping result of this 
study from a solid 4 x 4 to a solid 5 x 5 mesh is shown Table 5-1. The five avail- 
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able algorithms for solid target meshes show slight differences near the tran- 
sition zone of the two different fiber orientations of the original mesh. The 
element-based and the node-based results are slightly different since the fi- 
ber orientations information are element based with one integration point. 
Therefore, the node-based information is averaged. This averaging could lead 
to some inhomogeneities for some mapping algorithms. The Weighted Ele- 
ment formulation is especially developed for the usage in a CAE-chain with 
element based information [Mag16]. This approach is using the geometrical 
element information and weights them depending on the volume content. 
Due to this fact, the results are not continuous and therefore closer to the 
original mesh information. All other mapping algorithms are using a node 
based mapping, which leads to an initial deviation due to the averaging. 
Therefore, together with the other results of the study, the “Weighted Ele- 
ment” algorithm turned out to be the best. For further information regarding 
the different algorithms and their neighborhood search options it is referred 
to the MpCCl MapLib manual [Spi12]. 
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Table 5-1: Comparison of the different mapping algorithms from a solid 4x 4 to a solid 
5x5 mesh based on the first component (A11) of the fiber orien tation tensor 
Element Node 
torOcta 0 
E LT 
n 
= |: ”s 
wo as 
E= 
& = 
‘= 
je) o 
pom 
© 
> 
>] 
= 
3 
Ww 
o 
Q 
© 
Pet 
v 
om 
2 
~ 
© 
= 
3 
iL 
v 
Q 
© 
et 
v 
x 
= 
l meran 
Oo 
oc 
~ 
=o 
wo E 
o v 
= Lu 


134 


5.3. Linear elastic homogenization with anisotropic material cards 


5.3. Linear elastic homogenization with anisotropic 
material cards 


The import of the mapped result information and the homogenization is di- 
vided into two separate scripts to enable different homogenization methods 
within the Fibermap Plug-in. In the first step the import and the clustering are 
performed. This clustering prevents the generation of one material card for 
every single element. Due to this clustering, similar fiber orientations are con- 
flated in one group. To ensure the correct orientation of the global fiber ori- 
entation, a material orientation is assigned for every element. The clustering 
uses the eigenvalues (a1, a2, a3) of the fiber orientation tensor to find similar 
material properties, as the fiber orientation tensor can be expressed in its ei- 
genvector coordinate system as: 


Ay, Áz Aj3 ay 0 0 
= Anz Ass = a, 0 : (63) 
s a3 e1,82,83 


Due to the convention that the sum of the eigenvalues must be equal to one 
and the sequence of aı 2 a2 2 az, the following limits can be distinguished: 


0.33 0.5 1 
Aiso = (0.33 » Aplanar = (05) » Quni = (0) (64) 
0.33 0 0 


Due to these limits, the possible eigenvalue space can be reduced to the tri- 
angle plotted in Figure 5-7. The K-Means method is very well suited for group- 
ing point sets of this magnitude. The aim of this method is to solve the follow- 
ing optimization task: 


Min (Se din; € colli = Kell?) . (65) 


In this equation, C is the desired number of clusters, C, the set of points in 
the c-th cluster and the u. the midpoint of the c-th cluster. An algorithm to 
solve this problem is the so-called Lloyd's algorithm, which is implemented in 
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the Python library SciPy. Due to the Abaqus CAE implementation, this cannot 
be used within Abaqus. Therefore, a log file is used to solve this problem in 
an external program and to call the results from the list. 


© Eigenvalues 

O Isotropic Fibers 

© Unidirectional Fibers 
© Fiber in plane 

Œ Possible Positions 


3. principal axis 


1. principalaxis os 


1 05 


04 
0 02 va 


2. principal axis 


Figure 5-7: Clustering of material properties based on the eigenvalues of the fiber orienta- 
tion tensor. Nodes in the same clusters have the same color. [Mey16] 


All necessary input can be passed to the clustering script by using the GUI (see 
Figure 5-8). Within this GUI, the mapped result VTK file and the corresponding 
target part need to be selected. If a region of the target part should be omit- 
ted, the reserved element sets option could be used. This could be relevant, 
if the part consists of multiple Sections or part instances with different mate- 
rials. 


The number of the clusters can be defined in the dropdown menu of the “De- 
sired Number of Materials”. To ensure the accuracy as well as an acceptable 
computation time, 30, 100 or 200 clusters can be selected. If more than 200 
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clusters are used, the preprocessing is slowed down and the amount of re- 
quired memory during the simulation is increasing. 


| This Plugin imports mapped results (VTK) to Abaqus 


Mapped Result (*vtk): 


Choose Target Part: Model-1: zugprobe |v 


Log File(*csv): C:\Users\uadsy\Deskto| {Create File] 


Desired Number of Materials: |30 |r| 


Choose reserved element sets: | None 


Figure 5-8: GUI to import and cluster mapped results from a VTK file to Abaqus 


After the material clustering, the homogenization is done. Therefore, differ- 
ent homogenization approaches for stiffness and thermal properties are 
available, which are all based on linear-elastic material properties of the fiber 
and the matrix. 


For the stiffness homogenization, the Tandon-Weng approach is used to de- 
fine the fourth order Eshelby tensor S of unidirectional inclusions [Tan84]. 
This tensor describes the mechanical properties of a fiber-like spheroidal in- 
clusion under equal strain. Starting with the same properties for the matrix 
and the inclusion, a unidirectional two phase formulation is described by Tan- 
don and Weng [Tan84]. For the description of non-unidirectional composite 
materials, a tessellation over the hemisphere of possible fiber orientations is 
performed to generate an array of unidirectional anisotropic mechanical ma- 
terial properties. Within this approach, the fourth order fiber orientation ten- 
sor A is required. Therefore, the IBOF closure approximation (cf. Section 2.4 
and eq. 15) is used. To calculate the forth order anisotropic stiffness matrix C 
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based the tessellated unidirectional stiffness properties the fiber orientation 
tensor, the averaging approach from Advani and Tucker [Adv87] is used. 


The homogenization of the thermal properties follows the same process. For 
the generation of the unidirectional thermal properties, different methods 
depending on the used material system are available. The Laws approach 
[Law73] is based on a self-consistent theory and is developed for a multi- 
phase problem (n > 2) material and can describe anisotropic heterogeneous 
composites. When this approach is reduced to the two phases fiber and ma- 
trix, it is equal to the Rosen-Hashin method [Ros70], which is also available. 
Furthermore, the Scharpery method [Sch68] for mainly unidirectional and the 
Chamberlain [Cha68] approach for unidirectional transverse-isotropic com- 
posites is implemented. 


The linear elastic material and the thermal properties of the fiber and the ma- 
trix can be inserted to the homogenization via the GUI (see Figure 5-9). Fur- 
thermore, the scalar fiber aspect ratio and the fiber volume content have to 
be defined. Due to the object-based script programming and the prepared 
dropdown menu, further stiffness or strength homogenization approaches 
can be implemented easily. Based on the settings during the clustering, this 
script generates the chosen amount of anisotropic material cards and assigns 
the corresponding material orientation to every element. 
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This Plugin generates a material for each orientation cluster 


Clustering Result (*_labeled.csv): | C:\Users\uadsy\Deskto 


Choose Target Part: | Model-1: zugprobe kl 


Matrix Properties Fiber Properties 
Young's Modulus: | Young's Modulus: 
Poission Ratio: | Poission Ratio: 
Heat Expansion: | | Heat Expansion: 
Composite Properties 
Lenght/Diameter Ratio: 
Volume Fraction: 


Homogenization Scheme 


| Linear Elasticity: |Tandon Weng |y 


Expansion: | r 


Figure 5-9: GUI to define the material properties of matrix and fiber to homogenize the 
long fiber composite based on the fiber orientation tensor 


If results from structural simulations are used in further Abaqus simulations, 
the homogenized material properties can be exported into the VTK format. 
Therefore, the clustering result and the part instance can be selected in the 
export GUI (see Figure 5-10). This tool can also be used if Abaqus should not 
be used for the (thermo-) mechanical simulation. 
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This Plugin exports an Abaqus Part with homogenized properties to IRTG compliant VTK. 


Clustering Result (*_labeled.csv): | C:\Users\uadsy\Deskto 


Choose Part: |Model-1: zugprobe [x 


Figure 5-10: GUI to export homogenized material properties into the VTK format 


5.4. Effects ofthe consideration of fiber orientation 
in the structural simulation 


To show the effects of this CAE-chain, the results from the IRTG reference 
structure (cf. Section 4.3.2) are used. Therefore, two tension specimens are 
simulated, one with the mapped fiber orientation and the other with ahomo- 
geneous quasi-isotropic in-plane fiber orientation. The position of the tensile 
bar is given in Figure 5-11. Since the standard tensile specimen is slightly 
longer than the top surface of the beading, the clamping is considered iso- 
tropic where no mapped material is available. 
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Figure 5-11: Position of the tensile bar in the IRTG reference structure 


For the homogenization, 100 clusters are used. For both cases, the isotropic 
material properties given in Table 5-2 are used for the fibers and the resin. By 


assuming an averaged fiber diameter of 13 um and a fiber length of 25 mm, 


the aspect ratio is approximately 2000. The fiber volume content is set to 


40 %. 


Table 5-2: Material data for the homogenization 


Glass Fiber Resin 
Young’s Modulus [GPa] 73.0 3.1 
Poisson’s Ratio [-] 0.39 0.22 


By using the stiffness homogenization tool from Section 5.3, the material 


cards and the corresponding element sets are generated. To model the 1D 
load to the tensile bar specimen, the surfaces of the left clamping area are 
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constraint in all directions, while the right areas are only constrained in y- and 
z-direction. To apply the load in x-direction, a surface shear load is applied 
(see Figure 5-12). 


Figure 5-12: Boundary conditions and loads of the tensile bar specimen 


The results of the structural simulation are shown in Figure 5-13 with the von 
Mises stress distribution and in Figure 5-14 with the displacement. By com- 
paring the results from the mapped and the quasi-isotropic simulation, the 
mapped result shows slightly higher stress peaks and a slightly higher dis- 
placement. From this, it can be derived, that the average stiffness is lower 
compared to the quasi-isotropic sample. Furthermore, the mapped specimen 
shows non-symmetric local variations, resulting from the non-symmetric ge- 
ometry of the IRTG reference structure. This leads to anon-symmetric stress 
and deformation distribution. As these variations occur also over the thick- 
ness, a small out of plane torsion is predicted. 
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Figure 5-13: Von Mises stresses of the deformed tensile bars: mapped fiber orientation (top) 
and quasi-isotropic fiber orientation (bottom) from the IRTG reference 


structure 
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Figure 5-14: Displacement of the deformed tensile bars: mapped fiber orientation (top) and 
quasi-isotropic fiber orientation (bottom) from the IRTG reference structure 


This simple example shows the need of the consideration of manufacturing 
effects from the process simulation in the (thermo-) mechanical simulation. 
Furthermore, these effects need to be considered within the CAE-chain (see 
Figure 5-1) to optimize the process parameters (e.g. position of the initial 
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charge or the mold closing speed) or the geometry of the part within a topol- 
ogy optimization. 


5.5. Summary CAE Chain for long-fiber 
reinforced composites 


The process chain with the combination of the compression molding simula- 
tion and the structural simulation was shown for Sheet Molding Compound. 
Since also results from Moldflow can be used within the Fibermap Plugin, the 
presented CAE chain is also applicable to injection molding of thermoplastic 
and thermoset materials, as well as for the compression molding of long-fiber 
reinforced thermoplastics. 


By showing the different steps of the CAE chain with the export of results, 
converting to the neutral VTK-format, mapping, homogenization and struc- 
tural simulation, the user-friendly application was proven. The comparison of 
a homogeneous distributed fiber orientation and the mapped results from 
the process simulation at a tensile specimen show the need of such a CAE 
chain, since the maximal displacement differs by 7%. Since the consideration 
of the manufacturing induced fiber orientations leads to more realistic results 
in the structural simulation, global safety factors can be reduced, which leads 
to a reduced part weight and can be therefore the key to increase the usage 
of SMC. 
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For the mass production of fiber-reinforced polymers with Class-A surfaces, 
Sheet Molding Compound (SMC) is an outstanding material due to its good 
price to stiffness ratio and its very outstanding price to weight ratio. Due to 
the high possible degree of automatization and very low amount of waste due 
to a high material utilization factor, SMC is widely used in the automotive and 
the commercial vehicle industry. Nevertheless, new developments of new 
SMC formulations with higher fiber volume contents and sometimes with car- 
bon fibers extends the application area for SMC towards semi-structural parts 
[Hen17]. Due to these new applications, the holistic simulation of SMC parts 
becomes more important. Since commercial simulation tools are not suitable 
for the process simulation of SMC, new approaches need to be developed. 


Up to now, no standard is set to characterize the flow behavior of SMC to gain 
the material parameter for the process simulation. The latest developments 
were using simple compression rheometer with two pressure sensors to 
measure the pressure difference along the flow length [Dum07] [Gui12]. With 
this approach, the resin-rich lubrication layer at the mold surface and the rel- 
atively cold core region can be characterized under conditions close to the 
real process. The disadvantage of this available rheological tool is that multi- 
ple molding trials are necessary for one characterization run. Due to the in- 
homogeneity of SMC, this can lead to avoidable, method induced deviations. 
To overcome this disadvantage, a new rheological tool with a long flow path 
and seven pressure sensors was developed in the present work to obtain 
more reliable material data. Simultaneously, all press data are recorded. 


With this new rheological tool, four SMC materials were investigated. The first 
SMC material formulation is a typical Class-A SMC for exterior automotive 
parts. The second one is also a Class-A SMC formulation, where the replace- 
ment of a certain amount of fillers with micro hollow glass spheres led to a 
Low-Density Class-A SMC, developed by the Fraunhofer ICT. The third SMC 
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formulation is a semi-structural SMC, based on a Vinylester resin. This mate- 
rial has no fillers and an increased fiber volume content to carry higher loads 
in the final part. The fourth material is also a semi-structural SMC based on 
an Unsaturated Polyester Polyurethane Hybrid (UPPH) resin. This formulation 
contains no fillers and a further increased fiber volume content. To show the 
opportunities of the semi-structural SMC materials, this UPPH SMC was pro- 
duced with glass and carbon fibers. 


Due to the new design of the rheological tool and the characterization of dif- 
ferent SMC material formulations, the appearance of material compressibility 
could be detected along the flow length. Since Guiraud [Gui12] showed the 
compressibility under isothermal conditions for room-tempered SMC with an 
external lubrication film, the compressibility could now be shown under real 
process conditions. By quantifying this material behavior, the compressibility 
could be proven for the semi-structural SMC formulations over a large pro- 
cess window of closing speeds (constant velocity and constant nominal defor- 
mation rate). Contrary to this, the standard Class-A SMC formulations with 
fillers show an incompressible behavior. By using an empirical formulation, 
the compressibility was fitted over the process range to consider it in the ma- 
terial characterization models. Based on this compressibility model, the char- 
acterization model developed by Dumont et al. [Dum07] was extended to suit 
these semi-structural SMC materials. By using this new compressible model, 
which can also be used for incompressible material formulations, all SMC for- 
mulations except the glass fiber UPPH SMC were characterized. In a next step, 
the material parameters are compared with each other and with the data 
from the literature. Based on these data, some regularities could be derived. 
So it emerged that the power-law coefficients of the models for the lubrica- 
tion layer and the core region could be determined only depending on the 
material formulation. This reduces the characterization effort from four to 
two material parameters, which leads once again to a higher quality of the 
material data. Based on all this experience, this tool design with the corre- 
sponding compressible model has the potential to become a standard ap- 
proach to characterize SMC materials for advanced simulation models. 


146 


6. Conclusion 


Beside the material characterization for viscosity and compressibility, further 
new findings on material behavior could be made. Although the UPPH SMCs 
with glass and carbon fibers use the same resin, different flow behaviors were 
observed. While the carbon fiber UPPH SMC has a similar flow behavior like 
the other SMC formulations, the glass fiber UPPH SMC shows a different flow 
behavior. The standard SMC flow behavior is dominated by the thin resin-rich 
lubrication layer, which leads to an extensional deformation behavior in the 
core region. Since the pressure distribution of the glass fiber UPPH SMC with 
its pressure peaks at the flow front could not be described with such a flow 
behavior, optical evaluations and CT-scans were made to identify the reasons 
for the pressure peaks. Based on these evaluations, a probable flow model, 
which could describe this behavior, was developed. The basic idea of this 
model is that the lubrication layer consists of the outer sheet of the initial 
charge stack, which is flowing along the mold surface, while the core region 
is elongated only barely. This presented proposed flow model needs to be 
proven with further investigations. This is especially important since no fiber 
orientation could be predicted without a reliable flow model. This is one of 
the main issues when such a semi-structural SMC material is supposed to be 
used for automotive industry. 


The process simulation of SMC is the second focus of this thesis. Up to now, 
most process simulation models use a 2.5D shell approach based on a Hele- 
Shaw model with a varying thickness over the molding process. These models 
consider the lubrication layer and the core region by superimposing different 
models for both regions. Since this approach has problems when it comes to 
junctions, e.g. at ribs, these models are not suitable for nowadays needs. Fur- 
thermore, new developments in the area of combining SMC with local unidi- 
rectional fiber reinforcements require 3D approaches with a full Fluid-Struc- 
ture-Interaction. Therefore, commercial 3D SMC process simulation tools 
were compared with the experimental results. Even after a fitting of material 
parameter, these commercial tools were not able to predict the typical SMC 
pressure behavior with the high pressure increase at the flow front due to the 
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rheological stress which is followed by a threshold and the exponential pres- 
sure increase due to the hydrodynamic friction stress, while the material is 
flowing along the sensor. The differences between the experiments and the 
results from the commercial tool could be related to the used material mod- 
els. The commercial tools are using the Castro-Macosko [Cas82] model, which 
originates from the thermoplastic compression molding models and accounts 
for shear viscosity only. Therefore, neither the lubrication layer at the mold 
surface, nor the resulting deformation rate depending elongational behavior 
can be predicted. 


To overcome these disadvantages in the process simulation, a new approach 
based on the Abaqus Coupled-Eulerian-Lagrangian method was developed. 
This method offers a full Fluid-Structure-Interaction and can model large de- 
formation due to the combination of the Eulerian and Lagrangian approach. 
To enable this method to simulate the compression molding process of SMC, 
a user material subroutine (VUMAT) was developed, which contains the stress 
updates and the fiber orientation models. To calculate the current stresses in 
each increment, a superposition of the rheological stresses and the hydrody- 
namic friction stresses is performed. Therefore, the 2.5D models are ex- 
tended to full 3D ones. Additionally, a compression pressure due to the com- 
pression of the volume elements is considered. For the rheological stresses, 
a separate modelling of the extensional and shear viscosity is done, since the 
SMC needs to be considered as a Non-Newtonian compressible material. 
Therefore, the Trouton ratio [Tro1906] which couples the extensional and 
shear viscosity is not valid. For the extensional viscosity, a deformation rate 
depending power-law approach is extended by an Arrhenius function to con- 
sider the temperature dependent material behavior. For the shear viscosity, 
the state of art Castro-Macosko model is used. Since the lubrication layer is 
too thin to model it with a separate element layer within acceptable compu- 
tational time, it is modeled with a friction model. Therefore, a slip condition 
at the tool surface is used to ensure the correct velocity field of the core re- 
gion. Additionally, the hydrodynamic friction stresses are applied to the ele- 
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ments at the tool surface by using the velocity-dependent power-law ap- 
proach, suggested by Abrams et al. [Abr03]. To describe the fiber reorienta- 
tion during the compression molding simulation, the ARD-RSC model 
[Phe09a] is used. Within this model, a closure approximation to calculate the 
fourth order fiber orientation tensor out of the second order fiber orientation 
tensor is necessary. Therefore, the IBOF-closure [Chu02] is used, which is a 
compromise between accuracy and calculation time. Since not all necessary 
information are available within the VUMAT, several subroutines and scripts 
were written to pass all information from the initial configuration and the cur- 
rent time increment to the VUMAT. 


To validate this new process simulation model, the rheological tool from the 
characterization experiments was used. With the simulations of the different 
SMC material formulation, the typical SMC pressure pattern with the initial 
pressure increase and the following threshold could be predicted. Neverthe- 
less, these simulations revealed issues with the compressibility behavior of 
the SMC, which leads to errors in the simulation. Depending on the rheologi- 
cal stresses and due to the compressibility of the material, the material vol- 
ume is locally reduced. Due to numerical oscillations, this compression can 
localize dramatically, which leads to an extreme element deformation. Since 
the compression pressure model is not part of this work and further investi- 
gations need to be done, this issue was not solved. Therefore, a new model 
for the compression pressure needs to be developed, which might replace the 
scalar bulk modulus with a tensor modulus to consider deformation rate or 
velocity dependencies. Since this issues with the compression pressure are 
amplified within this model due to the reduced degrees of freedom (only U1 
and p), further tests with a 3D geometry were performed. Therefore, the ref- 
erence structure from the International Research Training Group was used. 
This reference structure also contains local unidirectional reinforcements, 
which were considered in the simulation model. The simulation results show 
the capability of the new SMC-specific process simulation approach. While 
the SMC is flowing in the cavity, the displacement of the local reinforcements 
could be predicted due to the full Fluid-Structure-Interaction within the CEL 
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method. For this simulation, the local reinforcements were considered iso- 
tropic and homogeneous. For further works, a more realistic unidirectional 
tape model should be used to predict the deformation more realistic, espe- 
cially if higher deformations are expected. 


To use the results from the process simulation for further simulations, the 
third focus of this thesis is to establish a virtual CAE-chain to transfer infor- 
mation between the process and the structural simulation. Therefore, the 
script-based Fibermap tool was developed, which can export finite element 
models and results of the process simulation to the neutral VTK-format, map 
information between different meshes and homogenize the material stiffness 
based on the fiber orientation tensors. The mapping is done by using the map- 
ping library MapLib, which is part of the MpCCl Multi-Physics Software Solu- 
tions provided by the Fraunhofer Institute for Algorithms and Scientific Com- 
puting (SCAI). MapLib is therefore integrated in the scripts for an easier usage. 
The stiffness homogenization step is divided into two sub-steps. In the first 
step, a clustering based on the fiber orientation tensors is done to avoid the 
creating of a material card for every element. By using the K-Means method, 
the possible fiber orientation tensors are clustered into 30, 100 or 200 groups, 
which makes the pre- and post-processing easier and reduces the initializa- 
tion time for the solver. Based on these clusters, which have one center fiber 
orientation tensor each, the homogenization is done. Therefore, a two-step 
scheme is used by combining Tandon-Weng approach [Tan84] for the unidi- 
rectional homogenization with the fiber tensor based averaging proposed by 
Advani and Tucker [Adv87]. Here, once again, the IBOF closure is used to cal- 
culate the fourth order Eshelby tensor for the ellipsoid inclusion of the fiber 
in the matrix material. Besides the stiffness homogenization, also a thermal 
expansion coefficient homogenization can be performed. To do so, the same 
two-step approach is used, but different homogenization schemes for the uni- 
directional homogenization are used. 


To show the effect of the CAE-chain, a comparison between the standard 
quasi-isotropic method and the mapped fiber orientation tensor dependent 
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material properties is performed. Therefore, the fiber orientation results of 
the IRTG reference structure are mapped onto a tensile bar. This simple ex- 
ample illustrates the differences due to a change of the elongation by 7%. By 
the consideration of the process simulation results with better material mod- 
els, the structural simulation comes closer to the experimental results and 
therefore, global safety factors can be reduced. This can be one key to in- 
crease the usage of SMC in the automotive industry. To further improve the 
CAE-chain in the context of the automotive demands, homogenization of 
strength and damage evolution need to be developed and added to the Fiber- 
map tool. 


Summarizing the three main focuses of this thesis, it could be shown that the 
reliable characterization, an SMC-specific process simulation model and the 
transfer of the results to the structural simulation, improves the prediction 
quality for the final SMC part. When neglecting one of these steps, the pre- 
diction quality will increase dramatically and part damages or too heavy parts 
are the result. Since both are knock-out criteria in the nowadays automotive 
development process, it is highly recommended to follow the presented pro- 
cess even while increasing the characterization effort and the computational 
time. 
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Greek Symbols 

Qi, 02 fitting parameter of relative volume change model 
Og gelation conversion 

B relative volume change 

ý magnitude of the strain rate tensor 

é strain rate tensor 

Ey volumetric strain rate 

n extensional viscosity 

Nps pure strain extensional viscosity 

K scalar reduction factor 

A hydrodynamic friction coefficient 

& ellipsoid particle shape function 

p density 

p time-related density derivative 

DER local thickness averaged axial stress 

xv thickness averaged in-plane viscous stress tensor 
o stress tensor 
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Nomenclature 


(033) volumetric averaged axial stress 

Ofric friction stress tensor 

Orheo rheological stress tensor 

Ov viscous stress tensor 

T* fitting parameter of the Castro-Macosko 


viscosity model 


@ mapping function from the initial material point 

$ mapping function from the reference Î to the 
spatial domain Q 

x motion of the material in the spatial domain Q 

X initial values of the particle motions in the 


reference domain 0 


y fiber orientation function 

Q spatial domain 

Qo initial configuration in the material domain 
a reference domain 


Latin Symbols 
A second order fiber orientation tensor 


A11, A22, A33 Fiber orientation in 1, 2 or 3 
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[Pa] 
[Pa] 
[Pa] 


[Pa] 


[Pa] 


[-] 


Latin Symbols 


directions respectively [-] 
Aint interaction area [m7] 
A fourth order fiber orientation tensor [-] 
B fitting parameter of the Castro-Macosko 
viscosity model [kg mt s] 
b width [m] 
b Left-Cauchy-Green tensor [-] 
bi-bs fitting parameters of the ARD model [-] 


Bı-Bßs fitting parameters of the IBOF closure approximation [-] 


c number of clusters [-] 
(0 ARD fiber interaction tensor [-] 
Ci, C2 fitting parameters of the Castro-Macosko 

viscosity model [-] 
D deformation rate tensor [s*] 
Do reference deformation rate [s+] 
Dc compressibility change deformation rate [s+] 
D” symmetric plastic parts of the velocity gradient [st] 
d diameter [m] 
F deformation gradient tensor [-] 
g penetration length [m] 
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Nomenclature 


ho 


hr 


lo 

def 
Imax 
lincomp 


Isx 


M 
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rheological function according to 
Dumont el al. [Dum07] 


closing speed 

height 

initial height of the IC stack 

finale plate height 

identity matrix 

bulk modulus 

fitting parameter of relative volume change model 
velocity gradient 


fourth order analytical functions of the 
eigenvalues of A 


length 

fitting parameter of relative volume change model 
initial length of the IC stack 

reference length of relative volume change model 
rheological tool length 

incompressible flow front position 


position of sensor x in the rheological tool 


fourth order analytical functions eigenvectors of A 


Latin Symbols 


Re 


To 


Tb 


Vo 


friction power-law coefficient 
rheological power-law coefficient 
surface normal vector 

thickness averaged pressure 

fiber orientation vector 

pressure 

body rotation matrix 

elastic part of the body rotation matrix 


distribution function of the components in a fourth 
order Tensor 


arbitrary reaction terms of the stress tensor 
reference temperature 


fitting parameter of the Castro-Macosko 
viscosity model 


time 

stretch tensor 
velocity 

Volume 

reference velocity 


vorticity tensor/spin matrix 


157 


Nomenclature 


Xo coordinates of the reference configuration 
XI, XII evaluating length in flow direction 

Xi coordinates of the current configuration 
Abbreviation 

ALE Arbitrary Lagrangian Eulerian 

ar aspect ratio 

ARD Anisotropic Rotary Diffusion 

CEL Coupled-Eulerian-Lagrangian 


D-SMC Direct-SMC 


EVF Eulerian volume fraction 
FAST Institute of Vehicle System Technology 
FLD field variable 


IAM-WK Institute for Applied Materials 


IBOF Invariant Based Optimal Fitting closure approximation 
IC initial charge 

ICT Fraunhofer Institute for Chemical Technology 

IRTG International Research Training Group 

LD Low density 
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[m] 
[m] 


[m] 


Abbreviation 


ORW3 


RSC 


SCAI 


SDF 


SDV 


SMC 


UD 


UP 


UPPH 


VE 


VUAMP 


VUEOS 


VUEXTERNAL 


VUFRICTION 


VUINTERACTION 


VUSDFLD 


VUMAT 


VUVISCOSITY 


wbk 


improved orthotropic fitted closure approximation for 
a wild range of interaction coefficient value 


Reduced Strain Closure 


Fraunhofer Institute for Algorithms and 
Scientific Computing 


Strain Reduction Factor 

state variables 

Sheet Molding Compound 

unidirectional 

Unsaturated Polyester 
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